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1.  Introduction 

Non-destructive  evaluation  of  remaining  life  and  reliability  of  structure  components  in 
case  of  aging,  cyclic  damage  and  cumulative  damage  is  critical  to  the  assessment  of 
reliability  of  the  entire  aging  aircraft.  Reliability  of  aging  components  can  be  estimated 
with  the  help  of  probabilistic  model  developed  on  the  basis  of  crack  propagator  concept 
CP  (probability  of  crack  formation  between  any  two  point)  and  Monte  Carlo  Simulation. 
Parameters  of  CP  can  be  evaluated  and  experimentally  verified  with  file  high  level  of 
precision  by  nonlinear  acoustic  technique  and  instrumentation.  This  instrumentation  uses 
the  interaction  of  ultrasonic  waves  with  vibration  (Nonlinear  Elastic  Wave  Spectroscopy 
-  NEWS)  to  collect  data  for  reliability  assessment  based  on  probabilistic  model  of  crack. 
In  Phase  I  principals  of  nonlinear  acoustic  approach  were  demonstrated  and  draft  design 
of  the  instrumentation  was  developed. 


The  overall  objective  of  the  Phase  I  proposal  was  to  develop  nonlinear  acoustic 
instrumentation  for  estimations  of  the  crack  parameters  and  to  use  date,  received  by  with 
this  instrumentation,  in  a  new  probabilistic  fracture  model  for  the  aging  aircraft  structural 
component.  The  specific  objectives  of  the  program  were  formulated  as  follows: 

1.  Propose  fracture  mechanism  model  for  the  materials  or  part  in  question. 

2.  Develop  nonlinear  acoustic  instrumentation  for  the  control  of  initiation  phase  of 
fatigue  damage.  Determination  of  correlation  between  nonlinear  acoustic  response 
and  crack  parameters. 

3.  Develop  the  theoretical  model  describing  the  modulation  of  acoustical  signal  by 
vibration  due  to  crack  and  relationship  between  the  level  of  modulation  and  the 
crack  parameters. 

4.  Investigate  the  application  of  the  “crack  propagator”  approach  for  estimation  of 

remain  time  life  of  aircraft  structural  component. 

5.  Verify  the  proposed  probabilistic  fracture  mechanics  model  for  the  material  or 
parts  in  questions. 

The  brief  results  of  phase  I  project  can  be  formulated  as  followed: 

1.1.  Probabilistic  approach  to  reliability  assessment  of  aging  aircraft  structural 
component  in  question 

The  probabilistic  approach  to  reliability  assessment  of  aging  aircraft  structural  component 
in  question  have  been  developed.  The  concept  of  Crack  Propagator  has  been  extended 
into  the  time  domain  by  incorporating  the  Crack  Layer  kinetic  theory  into  a  probabilistic 
framework  of  Statistical  Fracture  Mechanics.  This  model  allows  to  provide  numerical 
simulation  of  Crack  Layer  growth  from  an  initial  small  defect  to  the  ultimate  failure. 
Such  stimulation  provides  functional  relation  between  the  time  to  failure  and  the  initial 
defect  size  for  different  applied  load  and  structural  element  considered.  The  suggested 
model  is  used  for  calculation  of  conditional  reliability  function  that  provides  useful 
information  in  evaluation  the  life  expectance  of  structural  component  with  required  level 
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of  reliability.  Calculations  of  the  lifetime  expectance  as  a  function  of  normalized  applied 
stress  have  been  conducted. 


1.2.  Development  of  the  Nonlinear  Elastic  Wave  Spectroscopy  (NEWS) 

perimS  setup  to  observe  nonlinear  acoustic  interaction  between  ulu^und  and 
vibration  was  mounted.  The  instrumentation  provides  the  excitation  of  high-frequency 
ultrasonic  wave  and  the  low-frequency  vibration  in  the  material  under  mvestigahon  and 
measurement  of  the  interaction  product  -  modulation  of  sound  by  vibration, 
modulation  manifests  itself  as  a  sideband  component  in  the  spectrum  of  the  received 
ultrasound  signal.  The  instrumentation  consists  of: 

1  The  low-frequency  excitation  part  is  an  instrumental  hammer  with  attached 
piezoelectric  sensor  for  triggering  of  recording  signal.  The  several  mstmmen  a  s 
hammers  for  various  condition  of  vibration  excitation  were  build  and  tested. 

2.  The  high-frequency  ultrasound  part  mcludes  high  frequency  generator,  power  and 
piezoelectric  sensors  attached  to  the  investigated  part. 

3  The  registration  part  consist  of  piezoelectric  sensor  to  transfer  acoustical  signal  to 
electrical  signal  and  LeCroy  digital  oscilloscope.  The  pnmaiy  processing  was 
conducted  by  oscilloscope.  The  signal  was  also  transferred  to  computer  for 

secondary  processing. 

£  S"  ^  developed.  The  software  is  written  using 

MatLab  and  allows  to  conduct  the  following  operation  with  the  received  signals. 

a.  Spectral  analysis 

b.  Averaging  of  the  signal  due  to  impact  variation. 

c  Calculation  of  suggested  criteria  of  modulation  such  as  the  Nonlinear  Index  OSfl). 
The  modulation  criteria  in  the  frequency  domain  is  calculated  by  sununa  ion, 
normalization  and  averaging  of  the  modulation  spectra  for  the  certain  range  of 
ultrasonic  frequencies. 

1  4  NEWS  testing  of  different  materials  . 

The  preliminary  testing  was  conducted  for  different  kinds  of  materials  used  in  airsp 

industry.  The  following  materials  were  investigated: 

tSzrl1 used  in  F16  canopy  were  tested.  The  crack  growing  was 
Mtiated  under  cyclic  loading  .The  large  difference  in  the  modulation  signals  was 
observed  for  samples  with  small  cracks  in  comparison  with  samples  without  cracks  The 
dependence  of  the  nonlinear  index  on  the  crack  size  was  also  observed  and  investigated. 
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b.  Aluminum  strain  hardening  components 

The  aluminum  strain  hardening  components  (47x45x5  mm3)  were  tested  by  NEWS.  The 
modulations  were  observed  in  samples  with  different  crack  lengths.  The  nonlinear 
acoustic  signatures  for  different  samples  having  different  crack  size  were  obtained. 


c.  Adhesive  bounded  aluminum  plates 

The  nonlinear  acoustic  technique  was  applied  to  aluminum  bounded  plates  supplied  by 
Boeing.  The  samples  were  heat  treated  to  produce  some  kind  of  damage  and  reduce  the 
strength  of  bonded  joints.  The  experiments  had  shown  higher  level  of  modulation  for  the 
samples  having  the  low  level  of  strength  compare  to  normal  samples  that  were  not  acted 
by  treatment. 

d.  High  pressure  pipelines 

NEWS  was  applied  to  thin  high  pressure  pipelines  used  in  different  control  systems.  The 
damage  was  caused  by  multiple  tube  bending.  NEWS  registered  appearance  of 
nonlinearity  (accumulated  damage)  in  spite  of  no  visible  crack  was  observed. 

1.5.  Theoretical  model  describing  the  modulation  of  ultrasound  by  vibration  due  to 
cracks 

The  two  theoretical  model  describing  nonlinear  acoustic  effects  produced  by  a  single 
crack  as  well  as  a  cluster  of  microcracks  have  been  developed. 

a.  A  single  crack  model  considers  a  crack  as  a  domain  of  contact  between  two  rough 
surfaces.  The  applied  vibration  varies  surface  contact  area  within  the  domain  due  to 
deformation  of  asperities  in  contact,  leading  to  a  nonlinear  elastic  response.  The 
suggested  acoustical  crack  model  takes  into  account  the  resonance  properties  of  the 
investigated  specimens,  crack  size  and  crack  position.  The  expression  for  the  level  of 
side  components  in  a  spectra  of  the  simple  geometry  specimen  has  been  obtained.  The 
theoretical  model  predicts  the  amplitude  and  frequency  dependencies  of  the  side 
components  on  the  amplitude  levels  of  the  applied  signals  and  on  frequencies  of 
excitation.  The  linear  dependence  of  side  component  on  the  level  of  low  frequency 
vibration  was  determined.  This  dependence  allows  to  use  a  certain  modulation  criterion 
that  is  independent  on  the  level  of  excitation,  i.e.,  Nonlinearity  Index.  This  method  is 
especially  impotent  for  impact  produced  vibration  since  it  is  almost  impossible  to  provide 
the  constant  parameters  of  every  impact.  The  predicted  dependencies  of  side  components 
on  the  frequency  of  ultrasonic  waves  demonstrate  the  necessity  for  frequency  averaging 
of  the  side  components.  The  program  calculated  NEWS  output  data  for  parts  having 
simple  geometry  have  been  developed  and  the  model  calculations  have  been  conducted. 
They  demonstrated  the  dependencies  of  the  side  band  components  on  interaction 
frequencies  and  crack  position. 

b.  Model  of  multicrack  cluster.  This  model  considers  the  variation  of  the  nonlinear 
acoustic  parameters  in  the  material  with  microcracks.  The  suggested  approach  is  based 
on  Prof.  Chudnovsky's  work,  on  the  foundations  of  Statistics  of  Fracture  Statistic.  The 
developed  theory  connects  the  nonlinear  signature  of  the  material  with  the  crack 
distribution  characteristics  .  This  information  will  serve  as  a  basis  for  detection  of  the 
distributed  damage. 
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1.6.  Crack  location  by  NEWS. 

NEWS  pulse  location  method  is  based  on  the  radiation  and  reception  of  ultrasonic 
pulses.  The  difference  from  common  ultrasonic  pulses  scanning  method  is  an  application 
of  the  low  frequency  vibration  to  the  tested  component.  This  vibration  produces  the 
modulation  of  the  signal  reflected  from  cracks.  Experimental  verification  of  the  method 
has  been  done  on  steel  plates  having  different  kinds  of  defects  (a  hole,  an  open  crack,  a 
close  crack).  The  high  frequency  impulse  had  a  frequency  more  higher  than  1MHz,  and 
the  modulating  vibration  frequency  was  about  10  Hz.  It  was  demonstrated  that  this 
method  can  distinguish  a  crack  on  the  background  of  strong  reflection  from  other 
inhomogeneities  (open  holes  and  holes  with  bolts). 

1.7.  Mechanical  Design  of  the  Instrumentation 

Different  design  options  of  the  instrumentation  were  considered.  Kinematics  analysis 
was  performed  and  parameters  of  the  instrumental  hammer  were  defined.  Production 
drawings  of  the  instrumentation  are  under  consideration.  The  design  and  manufacture  of 
the  instrumentation  will  be  completed  in  Phase  II  of  this  project. 
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2.  Probabilistic  approach  to  reliability  assessment  of  aging 
aircraft  structural  component  in  question 

2.1.  Fundamentals  of  probabilistic  approach. 

The  primary  objective  of  our  efforts  is  the  development  of  a  computerized  system  for 
remaining  life  and  reliability  assessment  of  the  aircraft  structural  components  on  the 
basis  of  a  new  probabilistic  approach  combined  with  nonlinear  acoustic  characterization 

of  the  current  state. 

The  aircraft  aging  phenomena  usually  consists  of  three  stages: 

1  Submicro-  and  microdefect  accumulation  (incubation)  leading  to  crack  inaction. 

2.  Slow  crack  growth  (fatigue,  creep,  stress-corrosion  crack). 

3.  Dynamic  (unstable)  crack  propagation  leading  to  a  catastrophic  fracture. 

The  useful  lifetime  of  an  aircraft  structure  is  determined  by  the  first  two  stages,  since  the 
third  one  is  very  short.  The  formation  of  a  set  of  cracks  manifests  the  end  of  the  first 
stage  and  the  beginning  of  the  second  stage.  It  is  our  intention  to  utilize  a  nonlinear 
acoustic  technique  for  determination  of  the  beginning  of  the  second  stage  and  to  correlate 
the  basic  parameters  of  the  set  of  cracks  (the  distribution  of  crack  size,  number,  mutual 
position)  with  nonlinear  acoustic  response  of  the  structural  element. 

The  duration  of  the  second  stage  is  determined  by  fatigue  and  /or  corrosion  crack  growth 
rate  that  will  be  estimated  on  the  basis  of  Crack  Propagator  (CP)  concept.  CP  concept, 
recently  developed  within  Statistical  Fracture  Mechanics,  will  be  employed  for 
evaluation  of  the  remaining  lifetime.  The  Crack  Propagator,  i.e.,  the  probability  of  crack 
front  extension  from  one  position  to  another  during  a  specified  time  interval,  depends  on 
the  load,  structural  elements  geometry,  crack  size  and  fracture  surface  roughness. 
Reliability  of  aging  components  will  be  then  evaluated  with  the  help  of  probabilistic 
model  developed  on  the  basis  of  crack  propagator  concept.  To  fulfill  this  task  statistical 
distribution  and  failure  rates  will  be  experimentally  determined  for  different 
environmental  condition  of  crack  propagation.  Probability  density  function  for  the  stress 
and  strength  variation  over  the  time  will  be  evaluated  and  reliability  of  aging  components 
that  accounts  for  strength  and  toughness  deterioration  will  be  estimated. 

Statistical  characterization  of  the  current  crack  population  (size,  number  and  locations)  by 
nonlinear  acoustics  NDE  technique  is  one  of  the  challenges  of  the  proposed  method.  This 
technique  is  based  on  nonlinear  acoustic  principles  and  is  much  more  sensitive  to  fatigue 
damage  than  the  commonly  used  ultrasonic  NDE  technique.  The  foundation  for  the 
new  nonlinear  acoustic  NDE  method  is  the  significant  change  in  nonlinear  elastic 
response  of  the  structural  element  resulting  from  crack(s)  initiation  and  growth.  The 
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sensitivity  of  nonlinear  elastic  spectroscopy  (NEWS)  to  the  initiation  and  propagation  of 
cracks  is  orders  of  magnitude  higher  than  that  of  conventional  NDE  methods. 

The  nonlinear  acoustic  NDE  results  combined  with  the  information  regarding  geometry 
and  service  conditions  of  a  particular  structural  element  will  provide  a  sufficient  input 
data  for  evaluation  of  remaining  lifetime  and  reliability  by  Statistical  Fracture  Mechanics 
techniques. 

Basic  Concepts  and  Assumption 

Several  natural  assumptions  regarding  the  propagation  stage  of  fracture  process  underlie 
the  probabilistic  model: 

l  The  crack  path  is  randomly  selected  from  a  set  of  all  possible  (virtual)  crack  paths 
Q  (The  set  of  experimentally  observed  trajectories  depicted  in  Fig.  2.1  is  a 
pictorial  representation  of  a  sample  set  of  Q  ) .  Statistical  characterization  of  this 
set  is  based  on:  a)  fractographic  analysis  of  observed  crack  trajectories  for  the 
material  or  part  in  question,  b)  stress  analysis  and  c)  the  nonlinear  acoustic 
determination  of  the  size,  orientation  and  location  of  the  most  dangerous  defect 
(crack). 

2.  Crack  advance  along  a  particular  path  is  characterized  by  the  Crack  Layer  kinetic 
equations,  that  represent  a  mathematical  model  of  evolution  of  a  system 
consisting  of  a  crack  and  a  process  zone  (PZ),  formed  in  front  of  the  crack  tip. 

3.  Crack  advance  consists  of  a  sequence  of  elementary  steps,  controlled  by  Grifith 
criterion  Gi  >  2  y,  where  Gi  stands  for  the  energy  release  rate  (ERR)  due  to  crack 
extension  into  the  PZ,  and  y  =  y(0  is  the  specific  fracture  energy,  that  decreases 
with  time,  reflecting  material  degradation  within  the  process  zone. 

One  of  the  main  building  blocks  of  the  formalism  is  the  crack  propagator.  Let  us  consider 
a  component  under  a  load  containing  a  crack  coo  with  its  tip  as  x  (Fig.2.2).  First  we  define 
a  conditional  propagator,  P (x,t;X,T  /  co),  as  the  probability  that  the  crack  tip  advances 
from  position  x  at  the  time  t  through  position  X  at  the  time  T  >  t  along  the  path  o  ,  i.e.  the 
condition  Gi  >  2  y (t)  is  met  at  every  point  of  the  crack  trajectory  co. 

If  one  assumes  that  the  crack  trajectories  are  mutually  exclusive  (only  one  trajectory  is 
observed  in  every  individual  specimen,  then  the  probability  ,  P(x,t;X,T  )  of  crack 
extension  from  x  through  X  during  the  time  interval  T  -  t  can  be  written  as  total 

probability. 

P (x,t;X,T)  =  X  P (x,t;X,T /  (oK)  P(©k) 

Here  P(©k)  stands  for  the  probability  that  the  crack  “chooses”  a  path  coK  among  all  virtual 
paths  extended  from  x  to  X.  P(x,  t;X,  T )  is  called  the  crack  propagator  (CP).  In  other 


Fig.2.2.  A  sample  set  of  virtual  paths  leading  from  point  x  to  point  X 
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words,  CP  is  the  average  over  all  virtual  trajectories  of  the  conditional  probability  of 
crack  formation  between  the  points  x  and  X  along  a  particular  trajectory  ©K  during  the 
time  interval  T- 1. 

Once  the  CP  has  been  evaluated  for  a  particular  aircraft  structural  component,  the 
probabilities  of  crack  propagation  to  the  certain  depth,  crack  arrest,  distribution  of  critical 
load,  toughness  parameter,  time  to  failure,  and  reliability  function  can  be  readily 
computed. 

Evaluation  of  the  CP  results  from  three  major  tasks: 

1 .  Selection  of  a  set  of  virtual  crack  trajectories  Q  possessing  a  prescribed  roughness 
characterized  by  crack  diffusion  coefficient  and  the  fractal  dimension  of  crack 
trajectories. 

2.  Determination  of  the  conditional  propagator  P(x,  t;X,  T  /  coK) 

3 .  Computation  of  the  average  ?(x,  t; X,  T  /  ©K)  over  Q  by  means  of  a  functional 
integration. 

Monte-Carlo  stimulation 

A  numerical  approach,  most  convenient  for  development  of  a  software  for  reliability 
analysis,  is  Monte-Carlo  technique  [9].  It  allows  one  to  combine  the  physical  insight  and 
understanding  of  the  Statistical  Fracture  Mechanisms  with  the  flexibility  of  Monte-Carlo 
method.  In  particular,  it  allows  one  to  perform  parametric  studies,  necessary  for 
calibration  of  the  proposed  method. 

Following  the  spirit  of  Monte-Carlo  technique,  the  functional  integration,  i.e.  the 
averaging  over  the  set  of  virtual  crack  trajectories  Q  is  substituted  by  a  computer 
stimulation  of  a  large  number  of  virtual  crack  paths  followed  by  an  estimation  of  the 
frequencies  of  various  fracture  events  among  the  set  of  trajectories  Q.  The  approach  was 
successfully  applied  to  the  modeling  of  fracture  test  in  various  materials,  including 
different  composites  [10-14  ]. 

The  tasks  discusses  in  the  previous  section  can  be  reformulated  for  the  Monte-Carlo 
method  as  follows: 

1.  Computer  simulation  of  the  fracture  paths  with  a  prescribed  roughness,  i.e.,  the 
crack  diffusion  coefficient  and  fractal  dimension. 

2.  Numerical  solution  of  Crack  Layer  kinetic  equations  along  the  simulated 
trajectories  g>k  •  The  output  characteristics  for  each  trajectory  coK  from  the  set  of  { 
©K  }  are:  the  probability  distribution  of  the  depth  of  crack  penetration  during  time 
interval  T-t;  the  critical  crack  length;  the  time  t  *(coK)  of  crack  instability,  i.e., 
the  conditional  time  of  ultimate  failure  along  the  particular  crack  trajectory  <dks 
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conditional  reliability  R(  t  /  coK),  i.e.,  the  probability  of  “no  failure”  along  the 
particular  crack  trajectory  ooK  prior  the  time  t;  etc. 

3.  Statistical  averaging  of  die  output  characteristics.  Specifically,  the  reliability 
function  R(  t )  of  structural  component  in  question  is  defined  as  the  conditional 
reliability  R(  t  /  <dk  )  averaged  over  Q : 

R(  0  =  <  R(  *  /  ®k)  >  n 

for  the  Here  the  angle  brackets  <...>«  stand  for  the  procedure  of  statistical  averaging 
over  Q. 

Some  illustrative  examples  of  the  proposed  probabilistic  model  are  given  in  details  in 
[1,8,33]  and  will  serve  as  the  basis  for  the  reliability  assessment  particular  structural 
elements  under  specified  loading  and  environmental  conditions. 


2.2.  The  development  of  probabilistic  approach  for  reliability 
assessment. 

Kinetics  of  Fracture  Propagation. 

The  concept  of  Crack  Propagator  discussed  in  Phase  I  proposal  has  been  extended  into 
the  time  domain  by  incorporating  the  Crack  Layer  kinetic  theory  into  a  probabilistic 
framework  of  Statistical  Fracture  Mechanics.  The  fundamentals  of  Crack  Layer  theory 
and  it  implementation  in  probabilistic  setting  are  briefly  presented  below. 

The  crack  and  the  process  (damage)  zone  (PZ),  commonly  observed  in  front  of  the  crack  tip, 
are  together  considered  as  a  system  called  the  Crack  Layer  (CL)  [7-14,  33].  Since  the  size 
of  PZ  in  most  of  engineering  materials  is  much  smaller  than  the  crack  size,  the  CL  can  be 
characterized  by  two  parameters,  i.e.,  the  crack  and  PZ  sizes,  1  and  L  (Fig.2.3).  Following 
the  basic  principle  of  Thermodynamics  of  Irreversible  Processes,  the  driving  forces  of  the 
crack  and  PZ  growth  are  determined  as  derivative  of  Gibbs  potential  G  with  respect  to  the 
crack  and  CL  size  1  and  L  respectively: 
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Figure  2.3.  Illustration  and  schematics  of  the  Crack  Layer 


The  driving  forces  depend  on  properties  of  the  original  and  damaged  (within  PZ)  materials, 
and  account  for  the  energy  absorbed  within  the  PZ  during  its  formation,  as  well  as  on 
specimen  geometry  and  loading  conditions.  The  kinetics  of  CL  propagation  is  governed  by 
the  following  equations: 

d/dt l  =  ki  Jlcr,  d/dt  L  =  k2  Xpz  @.l) 


with  the  initial  conditions 

l  =  l0  and  L  =  L0  at  the  time  of  inspection  t  =  t0 


An  equilibrium  size  of  the  PZ,  Leq ,  is  determined  by  XPZ  =0. 

The  material  within  the  PZ  degrades  with  time  leading  to  a  decrease  of  specific  fracture 
energy  y  in  Grifith  condition,  i.e.,  y  =  y0  f(t),  where  yo  is  the  specific  fracture  energy  at  the 
time  to  and  f(t)  is  a  decreasing  function  of  time,  that  reflects  particular  mechanisms  of 
material  degradation  such  as  mechanical  fatigue,  electro-chemical  corrosion  or  thermal 
fluctuations,  to  name  a  few.  Kinetic  equations  (1)  together  with  the  PZ  material 
degradation  law  form  a  strongly  nonlinear  system  of  equations.  The  numerical  solution 
of  this  system  (the  fourth  order  Runge-Kutta  method)  reveals  two  basic  scenarios  of  CL 
behavior  shown  in  Figs.2.4.  The  CL  grows  continuously  (Fig.  6a),  when  during  the 
thermodynamic  forces  X,,  mdXL  are  positive  during  the  entire  process.  The  discontinuity  of 
the  CL  growth  (Fig.  6b)  appears  when  the  crack  driving  force  Xi  <  0,  so  that  the  crack  is 
arrested.  Crack  reinitiation  then  takes  certain  time  during  which  2}<r)  decreases  to  the  level 
of  Gi.  Than  the  crack  starts  to  propagate  through  PZ  until  it  reaches  an  undegraded 
material  becomes  arrested  and  stays  stationary  until  this  material  degrades  to  the  level  of  Gi 
to  allow  further  crack  advance.  Thus  the  next  cycle  begins.  Therefore,  a  discontinuous  CL 
growth  is  a  sequence  of  CL  stationary  states  and  relatively  fast  transitions  from  one  state  to 

another. 

Numerous  observations  of  various  authors  suggest  that  the  step-wise  discontinuous  CL 
growth  is  the  dominant  mode  of  crack  propagation  under  both  creep  and  fatigue  loading. 
The  average  rate  of  discontinuous  crack  growth  and  the  lifetime  are  mostly  determined  by 
the  duration  of  CL  stationary  configurations.  That,  in  turn,  depends  on  how  big  is  the  gap 
between  2y0  and  the  ERR  Gi.The  same  time,  the  lifetime  weakly  depends  on  the  values  of 
kinetic  coefficients  h  and  kL  that  affect  only  the  transition  time  from  one  stationary 
configuration  to  another.  This  time  interval  is  found  to  be  noticeably  shorter  than  the  CL 

stationary  states  (see  Fig.2.4b). 


15 


time,  f 


(a) 


(b) 


Figure  2.4.  Two  basic  scenarios  of  CL  growth 

Reliability  Function. 

Let  consider  an  acoustic  interrogation  of  the  structural  component  at  time  t0,  that  results 
in  detection  of  a  microcrack(s).  Apparently,  there  is  an  uncertainty  in  the  size,  location 
and  orientation  of  the  cracks.  Here,  for  an  illustrative  purpose  only,  we  assume  the 
location  and  orientation  of  the  crack  are  well  defined  and  only  the  size  is  uncertain.  A 
distribution  of  maximal  values  from  the  families  of  the  distribution  of  extremes  is  the 
distribution  of  choice,  since  the  largest  microcrack  is  most  likely  to  trigger  the  fracture 
process.  The  most  general  cumulative  distribution  function  F  (lo)  for  maximal  values  of 
initial  crack  size  lo  can  be  presented  as  follows  [8]: 
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The  range  of  possible  initial  crack  length  [/0min,/0max]  is  determine  from  the  acoustic  data 
and  the  scale  and  shape  parameters  a  and  (3  are  chosen  to  mach  as  close  as  possible  the 
histogram  of  the  experimentally  measured  crack  size. 

Numerical  simulation  of  CL  growth  in  single  edge  notch  specimen  from  an  initial  crack  l0 
to  the  ultimate  failure  provides  a  simple  functional  relation  between  the  time  to  failure  tf 
and  the  initial  crack  size  /0: 

tfU0)  =  A  +  B-  log(  /0) 

Here  the  coefficients  A  and  B  are  the  functions  of  applied  load  and  the  size  of  the 
structural  element  considered  [27].  Thus  the  cumulative  distribution  function  of  time  to 
failure  along  a  random  crack  trajectory  coK,  F(7f  /  coK  ),  results  from  strait  forward 
calculations  of  the  single  valued  function  t(l0)  of  the  random  variable  /0.  Therefore  the 
conditional  reliability  function  R(t/coK)is  also  determine,  since 

R(t/coK)  =  1  -  F(fr/<flK)  • 

Finally,  after  averaging  over  the  set  of  virtual  trajectories  Q  illustrated  on  Fig.2.1,  one 
finds  the  reliability  function  R(/)  =  <R(t/raK)>Q  for  various  level  of  applied  stresses 
depicted  on  Fig.  2.5.  The  calculations  have  been  performed  for  transparence  grade 
Polycarbonate  under  mechanical  stress  only  (aging,  temperature  variation  or  chemical 
aggression  have  not  been  considered).  However,  the  outlined  method  allows  to  account 
for  material  aging,  chemical  aggression  etc.  as  soon  as  the  kinetic  equations  for  the 
corresponding  processes  are  reported. 

The  reliability  function  thus  determined  provides  also  useful  information  in  evaluation 
the  life  expectance  of  structural  component  with  required  level  of  reliability.  This 
suggestion  is  illustrated  in  Figure  2.6  that  shows  the  lifetime  expectance  as  a  function  of 
normalized  applied  stress.  It  shows,  for  example,  the  lifetime  at  0.1  normalized  stress  is 
reduced  from  about  20  years  to  12  years,  i.e.,  by  40%  when  the  required  reliability 
increases  from  0.8  to  0.999. 

The  algorithms  of  these  types  of  computations  have  been  validated  on  various  examples. 


Figure  2.5.  The  reliability  function  for  various  level  of  normalized  applied  stress  ct/ct, 


Figure  2.6.  The  lifetime  dependence  on  the  normalized  stress  and  required  reliability  R 


18 


3.  Nonlinear  Acoustic  Technique  Principles 

There  are  several  nonlinear  effects  that  are  used  for  Nonlinear  Elastic  Wave  Modulation 
Spectroscopy  (NEWS),  i.e.  second  harmonic  generation,  modulation  of  sound  by 
vibration,  amplitude-dependent  internal  friction  and  resonance  frequency  shift  with 
amplitude. 

Historically,  the  first  applications  of  the  nonlinear  technique  were  performed  using  the 
second  harmonic  method.  This  method  involves  measurement  of  the  second  harmonic 
generated  by  the  nonlinear  distortion  of  a  sinusoidal  acoustic  signal  due  to  defects.  The 
first  test  of  this  method  in  the  United  States  was  made  as  early  as  in  1979  when  the 
second  harmonics  of  high  frequency  ultrasonic  waves  propagating  through  contact  was 
observed  by  Buck  et  al.  [2]  and  Richardson  et  al.  [3].  Another  modification  of  the 
second  harmonic  method  was  developed  for  the  quality  control  of  the  attachment  of 
thermal  protecting  tiles  for  the  Soviet  spaceship  "Buran",  Sutin  et  al.  [4].  One  more 
example  of  the  application  of  the  second  harmonic  method  is  the  detection  of  cracks  in 
large  (over  3  meter  long,  0.65  x  0.45  m2  cross-section)  graphite  electrodes  used  in 
aluminum  production  line,  Sutin  et  al.  [5]. 

The  team  from  Los  Alamos  National  Laboratory  has  developed  NDE  method  based  on 
the  measurements  of  the  resonance  frequency  shifts  and  damping  characteristics  as  a 
function  of  the  resonance  peak  acceleration  amplitude  [15,16].  This  method  is  called 
nonlinear  resonant  ultrasound  spectroscopy  (NRUS)  and  its  application  is  illustrated  in 
Figure  3.1. 
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Figure  3.1.  Experimental  procedure  for  nonlinear  resonant  ultrasound  spectroscopy. 
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A  rod  of  concrete  or  other  material  is  resonated  at  a  fixed  input  wave  volume  from 
frequency  a  to  frequency  b.  The  resonance  mode  lies  between  frequencies  a  and  b.  The 
input  wave  volume  is  increased  and  the  frequency  is  swept  from  a  to  b  again,  etc.  A  non- 
contact  laser  or  contact  accelerometer  is  used  to  detect  the  wave.  As  an  example,  an 
experiment  was  conducted  with  the  help  of  nonlinear  resonance  ultrasound  spectroscopy 
on  two  concrete  samples,  one  undamaged  and  one  damaged  by  the  alkali-silica  gel 
reaction  (ASR).  The  resonance  curves  show  a  significant  change  in  resonant  frequency 
as  drive  level  is  increased.  This  is  due  to  the  material  nonlinearity.  Other  nonlinear 
effects  including  harmonics  and  nonlinear  wave  dissipation  are  also  produced  (not 
illustrated).  Suffice  it  to  say  that  the  we  can  use  the  frequency  shift  to  characterize  the 
damage  in  the  material  using  the  “signature”  shown  in  the  figure. 


Undamaged  Damaged 


Figure  3.2.  Application  of  NRUS  to  differentiate  undamaged  and  ASR  damaged  samples 
of  concrete. 

The  illustrations  show  fundamental-mode  resonance  of  concrete  bars  for  seven 
progressively  increasing  volumes.  Note  the  very  different  nature  of  the  damaged  sample 
curves.  The  peaks  are  broader  and  the  shift  in  the  peak  with  volume  is  significantly 
higher  than  that  for  the  undamaged  sample.  The  peak  shift  and  the  peak  width  are  die 
nonlinear  diagnostics,  or  signatures. 

Another  approach  to  evaluate  the  nonlinear  acoustic  parameters  of  the  material  is  to 
measure  ultrasonic  waves  modulation  by  vibration,  [6,  21-28].  This  method  is  known  as 
nonlinear  wave  modulation  spectroscopy  (NWMS),  another  subset  of  NEWS.  The 
physical  nature  of  this  modulation  can  be  explained  with  a  simple  crack  model.  Consider 
the  sample  in  Figure  3.3  with  the  single  crack  shown  as  the  only  flaw.  The  method  uses 
the  interaction  of  a  high-frequency  ultrasonic  signal  with  low-frequency  vibration.  The 
low-frequency  vibration  signal  alternates  between  shrinkage  and  dilation  of  the  flaw. 
Consider  the  case  when  the  amplitude  of  the  vibration  during  the  compression  phase  is 
sufficient  to  close  the  crack  (Figure  3.3a),  and  the  subsequent  tension  opens  the  crack 
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(Figure  3.3b).  Suppose  the  ultrasonic  signal  has  a  frequency  at  least  an  order  of 
magnitude  larger  than  the  low-frequency  signal.  During  the  dilation  phase  of  the  low- 
frequency  cycle,  the  high-frequency  signal  is  partially  decoupled  by  the  opened  crack. 
This  reduces  the  amplitude  of  the  high-frequency  signal  passed  through  the  crack.  In  the 
other  half  of  the  low-frequency  cycle,  the  signal  produces  compression  and  the  crack  is 
closed.  The  closed  crack  does  not  decouple  the  ultrasonic  signal  and  the  amplitude  of  the 
transmitted  signal  increases.  This  results  in  an  amplitude  modulation  of  the  ultrasonic 
signal  as  shown  in  Figure  3.4;  i.e.,  low  amplitude  during  the  dilation  part  of  the  vibration 
signal  and  high  amplitude  during  the  compressive  phase.  Fourier  transformation  of  this 
signal  reveals  sideband  frequencies  that  are  the  sum  and  difference  of  the  frequencies  of 
the  ultrasonic  probe  and  vibration  signals.  These  new  frequency  components  are  the  sign 
that  a  flaw  or  crack  is  present. 


Figure  3.3.  Sample  with  flaw 

a.  closed  by  compression  b 

b.  open  under  extension 


Crack  opening 


Figure  3.4.  Amplitude  modulation  of 
probe  signal:  (a)  vibration, 

(b)  ultrasonic  signal 
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So  the  principle  of  the  NWMS  technique  based  on  ultrasound  modulation  by  vibration 
can  be  explained  by  using  the  following  scheme: 

Figure  3.5  shows  the  schematic  presentation  of  the  two  frequency  excitation  applied  to 
the  undamaged  part.  One  of  these  signal  has  frequency  ft  much  less  then  another  one  with 
frequency  f  2.  These  waves  passed  through  the  part  without  any  changes  in  frequencies. 

The  other  situation  is  observed  in  the  specimen  with  a  crack  (Fig.3.6).  The  low  frequency 
wave  modulates  the  high  frequency  wave  and  the  spectrum  of  modulated  wave  has 
additional  spectral  components.  The  spectrum  consists  of  primary  frequency  components 
fi  and  {2  and  from  two  additional  side  components  f+=fi+f  2  and  f_=fi-f2  .  These  side 
components  indicate  the  presence  of  crack. 


Figure3.5.  The  schematic  presentation  of  waves  passed  through  the  intact  specimen. 


Figure  3.6.  The  schematic  presentation  of  waves  applied  to  the  specimen  with  crack.  The 
interaction  of  waves  with  frequencies  fi  and  f2  and  leads  to  generation  of  side 
components  f+=fi+f2  and  f  =frf2  .  The  presence  of  these  side  components  indicates  the 

crack  presence. 
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4.  Experimental  Setup  and  Applied  Processing.  Polycarbonate 
Testing. 

As  it  was  described  in  previous  part  for  NEWS  we  have  to  apply  the  two  frequency  wave 
to  the  investigated  part  and  to  observe  the  side  components  due  to  modulation. 

There  are  several  ways  to  apply  two  frequency  waves;  one  is  to  use  two  separate  high  and 
low  frequency  transducers.  Simpler  way  is  to  use  one  ultrasonic  high  frequency 
transducer  and  to  produce  low  frequency  vibration  just  by  impact.  Tide  impact  by 
hammer  excites  the  free  vibration  of  inverstigated  part  and  these  vibration  gives 
modulation  of  probe  ultrasonic  wave. 

The  experimental  setup  to  observe  nonlinear  acoustic  interaction  between  ultrasound  and 
vibration  has  to  provide  the  high-frequency  ultrasound  excitation  wave  and  the  low- 
frequency  vibration  in  the  material  under  investigation,  and  measure  the  interaction 
product  —  modulation  of  sound  by  vibration.  This  modulation  manifests  itself  as  a 
sideband  component  in  the  spectrum  of  the  received  ultrasound  signal.  The  scheme  of 
typical  experimental  setup  is  presented  in  fig.  4.1.  The  instrumentation  consists  of: 

1.  The  low-frequency  excitation  part  -  instrumental  hammer  with  attached 
piezoelectric  sensor  for  triggering  of  recording  signal. 

2.  The  high-frequency  ultrasound  part  includes  high  frequency  generator,  power 
amplifier,  and  piezoelectric  sensors  attached  to  the  investigated  part. 

3.  The  registration  part  consist  of  piezoelectric  sensor  to  transfer  acoustical  signal  to 
electrical  signal  and  LeCroy  digital  oscilloscope  (Waverunner  LT322),  The 
primary  processing  was  conducted  by  oscilloscope.  The  signal  was  also 
transferred  to  computer  for  secondary  processing. 


Fig.4.2.  presents  the  photograph  of  the  actual  experimental  setup. 
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Figure  4.1.  The  experimental  setup  for  NWMS  measurements. 

For  a  nonlinear  acoustic  test  the  samples  were  supported  by  soft  foam  or  were  hanged  by 
thin  threads.  Two  ultrasonic  transducers  were  glued  to  the  surface  of  each  sample;  one 
was  used  as  a  transmitter  and  another  as  a  receiver  of  high  frequency  ultrasound  wave. 
Low  frequency  vibration  was  generated  with  an  instrumental  hammer,  which  produced 
calibrated  and  repeatable  impacts  of  controllable  amplitudes.  The  hammer  was  equipped 
with  an  accelerometer  used  for  synchronization  with  the  data  acquisition  system. 

The  signal  from  accelerometer  attached  to  the  hammer  is  shown  in  Figure  4.3.  The 
sensitivity  of  accelerometer  with  preamplifier  was  0.01  mV/  m/sec2,  it  means  that 
maximum  of  acceleration  during  impact  reaches  33000  m/sec2.  Taking  into  account  that 
the  hammer  mass  is  93  g  the  force  produced  by  hammer  is  about  3000  N. 


Figure  4.2.  Photograph  of  actual  experimental  setup. 


Voltage,  V 


Figure  4.3.  The  snapshot  of  impulse  produced  by  hammer. 

The  effect  of  ultrasound  modulation  by  vibration  can  be  clear  observed  in  time  domain 
for  spaciman  having  relatively  large  size  of  defects.  One  of  such  parts  is  Automoble 
Bearing  Cap  presented  in  fig.  4.4.  One  of  bearing  caps  has  a  crack  with  the  highest  level 
of  nonlinearity  among  all  samples.  The  time  series  of  the  received  signal  after  impact  for 
high  frequency  excitation  205  kHz  is  shown  in  Fig.  4.5.  It  can  be  seen  that  the 
modulation  takes  place  just  after  the  impact. 


Figure  4.4.  The  picture  of  Automobile  Bearing  Cap 


Figure  4.5.  The  time  series  of  the  detected  signal  signal  from  the  brearing  cap  with 
crack  (  frequency  205  kHz).  The  moment  of  impact  is  shown  by  arrow  and  it  is  seen 
the  amplitute  modulation  of  the  received  sinal  after  impact. 


Spectral  amplitude 


The  two  spectra  obtained  by  averaging  of  10  impacts  for  the  intact  and  cracked  bearing 
cap  and  the  bearing  cap  are  presented  on  the  fig.4.6  and  fig.4.7.  The  frequency  of 
excitation  was  192  kHz.  A  similar  difference  in  spectra  is  seen  for  any  frequency  studied. 
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Fig.4.6.  The  spectrum  of  the  received  signal  after  impact  for  cracked  bearing  cap.  There 
are  many  side  components  due  to  ultrasound  modulation  by  vibration. 
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Fig.4.7.  The  spectrum  of  the  received  signal  after  impact  for  intact  bearing  cap.  There  is 
only  main  frequency  component  and  no  side  components. 


The  results  of  modulation  depends  on  the  applied  ultrasonic  frequency,  sensor  position. 
The  way  to  avoid  such  influence  is  summation  average  proceducer  applied  to  the  the 
signal  with  varied  frequency. 

We  will  illustrate  the  summation  average  procedure  on  the  example  of  NWMS  technique 
application  for  crack  detection  is  polycarbonate  sample. 

The  specimens  were  machined  from  the  polycarbonate  used  for  aircraft  canopies.  The 
specimen  geometry  is  shown  in  Figure  4.8.  The  experimental  setup  presented  in  fig  4.1 
and  4.2.  was  used.  The  fatigue  cracks  were  developed  during  four-point  bending  cycling 
regime.  The  crack  lengths  are  given  in  Table  4. 1 . 


Table  4.1.  Tests  results.  Polycarbonate  specimens;  0.5  inches  thickness;  Four-point 
bending  cycling  regime. 


Specimen  Number 

Crack  length,  mm 

Number  of  cycles 

Stress,  psi 

1 

Reference  specimen 

0 

0 

2 

2 

85000 

10000 

3 

6.3 

128000 

10000 

4 

Initiation 

32000 

10000 

Figure  4.8.  Polycarbonate  specimen  with  attached  sensors. 


The  spectra  after  impact  for  different  samples  and  frequencies  are  presented  in  fig.  4.9.- 
4.15.  The  fig.4.9.  shows  the  typical  spectrum  for  undamaged  specimen  having  low  level 
of  side  components.  Figures  4.10  -  4.13  present  spectra  for  cracked  specimen  for 
different  frequencies  of  ultrasonic  wave.  It  is  seen  that  the  received  spectra  moved  along 
frequency  axis  and  these  spectra  can  not  be  averaged  directly.  We  have  developed  the 
summation  averaged  procedure  that  allows  to  average  side  components  for  different 
ultrasonic  frequencies.  This  procedure  is  working  by  the  following  way: 

Once  the  impact  device  strikes  the  object,  a  high-pass-filtered  time  signal  is  captured  on 
the  digitizing  device.  This  signal  contains  the  information  required  to  discern  whether  or 
not  damage  exists  in  the  test  sample.  In  the  simplest  case,  a  Fourier  transform  is  taken  of 
the  time  signal,  and  modulation  sidebands  are  looked  for  around  the  high  frequency 
ultrasonic  signal.  This  signal  is  stored  on  the  computer  for  further  analysis.  A  nonlinear 
effect  from  crack(s)  in  the  test  object  creates  modulation  of  ultrasound  by- vibration  and 
side  band  component  present  in  the  sample.  The  levels  of  sidebands  in  undamaged 
objects  are  much  less  then  in  damaged  object. 
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The  next  step  in  the  experimental  procedure  is  to  change  the  pure  tone  frequency(ies), 
and  repeat  the  experiment.  This  is  called  a  step-frequency  sweep.  This  procedure  is 
reneated  for  some  number  of  different  pure  tone  frequencies,  one  at  a  time.  For  instance, 
one  could  use  eleven  pure  tones  successively  at  185kHz,  186kHz,  187kHz....  195kHz. 
(The  spectra  for  frequencies  185,  190,  192  and  195  kHz  are  shown  in  fig.  4.10-4.13).  At 
each  pure  tone  frequency,  the  impact  device  is  applied  and  the  detected  signal  is 
collected.  The  spectra  of  any  signal  are  transformed  to  low  frequency  band  such  matter 
that  main  frequency  is  transformed  to  zero.  After  such  transformation  a  spectrum  of  any 
signal  will  consists  of  main  component  having  frequency  zero  and  side  components  with 
positive  and  negative  frequencies  around  main  frequency.  These  transformed  spectra  can 
be  averaged  and  the  result  does  not  depend  on  ultrasonic  high  frequency  signal. 

The  results  of  the  such  averaging  are  presented  in  Figures  4.14-4.17  for  all  four  tested 
samples.  The  averaged  level  of  the  side  components  is  the  higher  for  the  sample  with 
biggest  crack  (fig.4.15)  and  becomes  less  to  the  sample  with  smaller  crack  (fig.4.16).  The 
level  of  side  components  for  polycarbonate  without  crack  is  much  less  (  see  fig.  4.14  and 

4.17). 


Figure  4.9.  The  spectrum  of  the  ultrasonic  signal  with  frequency 
190  kHz  for  sample  #0  (undamaged). 


Figure  4.10  The  spectrum  of  the  ultrasonic  signal  with  frequency 
185  kHz  for  sample  #1  (6.3  mm  crack). 
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Figure  4.12.  The  spectrum  of  the  ultrasonic  signal  with  frequency 
192  kHz  for  sample  #1  (6.3  mm  crack). 
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Figure4.13.  The  spectrum  of  the  ultrasonic  signal  with  frequency 
195kHz  for  sample  #1  (6.3  mm  crack). 
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Figure  4.14.  The  spectrum  of  the  normalized  and  averaged  signal 
for  sample  #0  (undamaged). 
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Figure4.15.  The  spectrum  of  the  normalized  and  averaged  signal  for 
sample  #1  (6.3  mm  crack). 


35 


Frequency,  Hz 


Figure  4. 1 6  The  spectrum  of  the  normalized  and  averaged  signal 
for  sample  #2  (2mm  crack). 
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Figure  4.17.  The  spectrum  of  the  normalized  and  averaged  signal 
for  sample  #3  (undamaged). 
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This  summation  average  procedure  is  used  to  improve  the  signal  to  noise  ratio  and  to 
inspect  the  results  of  many  experiments  simultaneously.  The  result  of  the  summation 
average  provides  a  yes/no  result.  Is  the  object  damaged  or  not?  The  number  of 
frequencies  and  the  number  of  steps  in  frequency  are  not  critical.  The  more  frequency 
steps  the  better  the  signal  to  noise  ratio  is  in  the  average  of  the  Fourier  spectrum,  and  the 
more  robust  the  measurement  is  because  more  and  more  sidebands  are  used  in  the 
decision.  Because  die  averaging  is  done  in  the  frequency  domain,  the  phase  of  each 
successive  time  signal  collected  in  unimportant. 
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5.  Results  of  Experimental  Investigations 

This  section  of  report  presents  applications  of  NEWS  technique. 

5.1.  Observation  of  the  modulation  of  ultrasonic  waves  by  impact  produced 
vibration  in  strain  hardening  components 

This  example  demonstrates  how  the  developed  technique  works  for  strain  hardening 
aluminum  components.  These  components  have  size  47x45  x5  mm3  and  the  picture  of 
component  is  shown  in  fig.  5.1. 


Figure5.1.  Picture  of  the  stress  hardening  component. 

The  crack  was  grown  by  cyclic  loading.  The  crack  had  arisen  in  the  tip  of  notch  and  this 
crack  was  hardly  seen;  we  had  to  use  microscope  to  observe  the  crack.  The  length  of 
crack  was  from  1.5  to  3  mm. 

Four  damaged  and  four  intact  samples  were  tested.  The  cracked  samples  were  marked  as: 
#  55,  #  47,  #  65,  #  43  and  the  intact  sample  were  marked  as  Al,  A2,  A6,  A7. 

The  experiment  was  carried  out  with  small  PZT  piston  sensors  glued  to  the  flat  surface  of 
samples(  See  fig.5.1  and  5.2.).  The  samples  were  hanged  by  a  fine  thread  (wire). 

The  experimental  setup  is  presented  in  fig.5 .2. 
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Figure5.2.  The  experimental  setup  to  investigate  stress  hardening  components. 
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As  it  was  described  above  the  presence  of  crack  is  indicated  by  interaction  between  high 
frequency  ultrasonic  wave  and  low  frequency  vibration  produced  by  hammer  impact. 

The  sample  vibration  produced  by  impact  was  picked  up  by  the  same  piezoelectric  sensor 
and  were  recorded  by  oscilloscope.  The  snapshot  of  low  frequency  vibration  signal  after 
impact  is  presented  in  fig.  5.3  and  the  spectrum  of  this  signal  is  shown  in  fig.  5.4. 


Figure5.3.  The  time  sequence  of  vibration  produced  by  impact 


Figure5.4.  The  spectrum  of  impulse  produced  vibration 
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Two  waveforms  of  ultrasonic  signals  after  impact  for  frequency  203  kHz  for  cracked  and 
uncracked  samples  are  given  in  fig.5.5  and  5.6.  We  can  observe  modulation  of  ultrasonic 
wave  after  impact  in  the  cracked  sample  and  there  is  no  modulation  for  the  sample 
without  crack  (see  fig.5.6). 

The  more  quantitative  analysis  have  been  conducted  by  using  the  summation  averaged 
procedure  described  in  previous  part.  This  procedure  was  applied  to  the  signal  recorded 
in  the  frequency  band  200-210  kHz  for  set  from  8  samples  and  the  results  are  presented  in 
fig.  5.7-5.14. 

The  large  difference  of  averaged  spectra  for  damaged  and  undamaged  samples  can  be 
clearly  seen. 


Fig  5.8.  Averaged  summation  spectrum  for  sample  #43,  cracked 
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Figure  5.9.  Averaged  summation  spectrum  for  sample  #A2,  uncracked 
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Figure  5.10.  Avareged  summation  spectrum  for  sample  #47,  cracked 
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Figure  5.12.  Averaged  summation  spectrum  for  sample  #55,  cracked 
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Figure  5.13.  Averaged  summation  spectrum  for  sample  #A7,  uncracked 


Figure  5.14.  Averaged  summation  spectrum  for  sample  #65,  cracked 
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The  summation  averaged  procedure  was  applied  to  the  signal  recorded  in  the  frequency 
band  190-196  kHz  and  the  results  of  are  presented  in  fig.  5.17-5.18. 

It  is  seen  side  components  with  frequency  about  800  Hz  in  the  damaged  sample  and  there 
are  not  side  components  in  the  intact  sample.  The  difference  in  the  modulation  is  not  so 
clear  as  in  previous  tests 


Figure  5.17.  Summation  averaged  spectrum  for  the  damaged  sample  #4 


Figure  5.18.  Summation  averaged  spectrum  for  the  undamaged  sample  #  1 
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5.3.  Application  of  NEWS  for  quality  assessment  of  bending  tubes 

The  test  for  bending  tubes  demonstrate  that  the  developed  technique  can  be  applied  not 
only  for  crack  detection  but  also  for  quality  assessment  when  crack  did  not  exist. 

The  similar  technique  was  applied  for  thin  steel  tubes  used  in  aircraft  industry.  The 
picture  of  tube  is  presented  in  fig.  5.19. 


Figure  5.19.  Picture  of  tube 

To  produce  damage  some  of  tubes  were  affected  by  bending  up  to  90  degrees  and  back. 
The  numeration  of  tubes  and  parameters  of  bending  are  presented  in  table  5.1. 
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Table  5.1.  Parameters  of  tested  tubes. 


Numeration 

Length,  mm 

Number  of 
Bending 

A1 

170 

10 

A2  n 

165 

1 

A3 

170 

5 

B1 

210 

3 

B2 

200 

1 

B3 

200 

25 

B4 

205 

6 

B5 

215  1 

20 

B6 

205  ~~ 

0 

Cl 

138 

0 

C2 

138 

25 

The  results  of  summation  averaging  are  shown  in  figures  5.20  -. 

It  is  seen  that  the  tubes  bended  20  and  more  times  (#B3,  B5  and  C2)  have  distinguish 
modulation  (  clear  side  band  components(  while  other  tubes  have  no  modulation. 

The  observed  results  are  the  first  experimental  verification  to  for  damage  detection  in 
small  tubes  and  this  method  can  be  quickly  improved  be  choosing  the  most  sensitive 
frequency  range  and  conditions  for  impact  low  frequency  excitation 


Figure  5.21.  Summation  averaged  spectrum  for  tube  #  A2  ( 1  bending) 


Figure  5.24.  Summation  averaged  spectrum  for  tube  #  B3  (  25  bending) 
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Figure  5.25.  Summation  averaged  spectrum  for  tube  #  B4  ( 6  bendings) 


Figure  5.26.  Summation  averaged  spectrum  for  tube  #  B5  ( 20  bendings) 
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Figure  5.27.  Summation  averaged  spectrum  for  tube  #  B6  (  0  bending) 
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Figure  5.28.  Summation  averaged  spectrum  for  tube  #  Cl  (  0  bending) 
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Figure  5.29.  Summation  averaged  spectrum  for  tube  #  C2  (  25  bendings) 


55 


6.  Theoretical  Models  Describing  the  Modulation  of 
Ultrasound  by  Vibration  due  to  Defects 

6. 1.  The  single  crack  model 

The  developed  Vibro-Acoustic  Modulation  NDE  technique  deals  with  contact-type 
defects,  such  as  cracks,  unbonds,  delaminations,  etc.  The  physical  nature  of  the  defect- 
related  nonlinearity  can  be  most  easily  understood  when  we  consider  two  simple  models 
of  contact-type  defects.  The  first  model  was  suggested  by  Richardson  [3].  He  considered 
a  crack  or  contact  defect  as  planar  interface  separating  two  semi-infinite  elastic  materials, 
with  the  surfaces  in  intimate  contact  but  having  no  traction  forces  across  the  interface.  As 
a  longitudinally  polarized  elastic  wave  applied  to  the  interface,  the  surfaces  will  stay 
together  and  move  as  one  during  the  compression  phase  of  the  wave  and  then  separate 
under  the  tensile  phase.  Such  a  periodic  separation  causes  distortion  of  the  passing  wave 
(nonsymmetry  of  the  compression  and  tensile  phase)  and  leads  leading  to  the  high 
harmonics  formation. 

The  second,  more  realistic  model  of  a  defect  is  a  contact  between  two  rough  elastic 
surfaces  [29,  30].  The  applied  stress  will  vary  the  contact  area  within  the  defect  due  to 
deformation  of  asperites  in  contact,  leading  to  a  nonlinear  elasticity  of  the  defect.  Such  a 
contact  behaves  as  a  nonlinear  spring  whose  stiffness  is  proportional  to  the  contact  area 
within  the  defect  interface.  The  linear  model  of  the  spring-type  crack  interface  [31] 
connects  displacements  and  stresses  on  both  sides  of  an  interface  using  spring  stiffness 
constant  K. 

According  to  this  model  the  internal  stress  Acr  is  normal  to  the  crack  interface  is  given  as 
function  of  crack  deformation  by  (6.1): 


Acr=K%,  (6.1) 

where  £  is  the  variation  of  the  crack  thickness,  £  =  U+  -  U. ,  U±  are  the  deformations  of 
the  opposite  sides  of  the  crack.  This  simple  model  can  be  extended  to  include  nonlinear 
behavior  of  the  interface  introducing  dependence  of  the  spring  stiffness  on  the  crack 
thickness,  K=  K(Q.  It  means  that  the  applied  vibration  effectively  changes  the  stiffness 
of  the  crack. 

Using  the  Taylor’s  expansion  of  the  spring  stiffness  dependence  and  leaving  only  the 
first-order  term,  since  stresses  in  acoustical  signals  used  in  NDT  are  rather  small,  we 
obtain: 

K-Ko+aZ, 


Aa  =  (Kq+  a%)  g , 


and 


(6.2) 
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where  a  is  the  nonlinear  coefficient  (a£;  «K0).  This  nonlinear  model  of  a  defect  allows 
for  description  of  various  nonlinear  effects,  such  as  generation  of  second  harmonic  and 
modulation  of  ultrasound  by  vibration.  The  nonlinear  coefficient,  a,  can  be  used  for  the 
quantitative  description  of  these  effects. 

Some  aspects  of  the  sound  modulation  by  vibration  can  be  understood  using  theoretical 
model  of  a  plane  crack  in  a  thin  bar,  as  shown  in  Fig.6.1.  As  a  nonlinear  element,  the 
defect  can  transform  part  of  the  incident  acoustic  energy  into  acoustic  waves  with 
different  frequencies  (harmonics  and  combination  frequencies),  effectively  becoming  a 
source  of  the  nonlinear  acoustic  waves.  In  the  case  of  modulation  of  ultrasound  by 
vibration,  we  consider  a  flaw  as  a  simple  plane  source  of  the  acoustic  waves  with 
combination  frequencies  (o3=  atj  ±co2,  where  is  the  frequency  of  applied  vibration  and  co2 
is  the  frequency  of  the  incident  ultrasonic  wave. 
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Figure  6.1 .  The  considered  model  of  the  crack  in  solid  bar 

For  simplicity,  we  analyze  interaction  of  waves  propagating  only  in  x-direction  having 
wavelength  of  vibration  greater  than  bar  thickness  H«  L,  where  L  is  the  bar  length,  so 
die  thin  bar  approximation  can  be  used.  The  crack  interface  is  normal  to  the  direction  of 
wave  propagation. 

Let  us  consider  the  interaction  at  the  crack  interface  of  low  frequency  flexural  vibration 
(modulating  vibration)  having  longitudinal  deformation 

uj  (x,z,t)  =  U](x,z)  cos  (o\  t  (6.3) 

with  high  frequency  longitudinal  ultrasonic  which  having  deformation 

U2  (x,z,t)  =  U2(x,z)  cos  a>2 1  (6.4) 

We  assume  that  the  crack  has  an  effective  thickness  D  and  length  h  which  are  smaller 
than  the  wavelengths  of  the  incident  waves,  so  the  crack  does  not  have  significant  effect 
on  the  distribution  of  strain  and  stress  fields  in  both  low  and  high  frequency  waves. 


The  applied  stresses  change  the  thickness  of  the  crack,  so  the  crack  thickness  variation 
under  the  applied  stresses  can  be  determined  as 


r 


£=  £)  =  £)  EEj _cos  (0^  +  D 

dx  dx 


dU  2 

- —  COS  CD  9 1  - 

dx  2 


(6.5) 


The  nonlinearity  of  the  crack  is  characterized  by  the  coefficient  a  defined  by  (6.2).  This 
nonlinearity  leads  to  interaction  (modulation)  of  the  higher  frequency  ultrasonic  wave  by 
the  lower  frequency  vibration.  As  a  result  of  such  modulation  the  sidebands  or  wave  with 
combination  frequencies  will  be  formed  (0)3=  coj  ± 0)2 > 


Substitution  of.(6.5)  into  6.2)  gives  an  expression  for  the  stress  Acr3  at  the  combination 
frequencies: 


4  „  dU  r(xnzn)  dU ?( xnz0 ) 

A g3=  ccd2  — 1  -  — 21  0>  0  cos  03  t , 

dx  dx 


(6.6) 


where  x0  z0  are  the  coordinates  of  the  crack.  This  stress,  applied  to  the  surrounding  crack 
material  produces  the  acoustical  wave  with  sum  and  difference  (combination) 
frequencies. 


Let  us  consider  standing  wave  flexural  vibrations  (resonance  vibration),  which  are 
common  for  many  practical  applications.  This  situation  also  takes  place  in  the  experiment 
described  below. 


The  standing  low  frequency  flexural  waves  in  thin  elastic  bar  are  determined  by  the 
formula : 

w=  A  coskj  xcoscojt .  (6.7) 

Here  w  is  the  deformation  of  the  bar  in  z  direction,  A  is  the  amplitude  of  the  deformation. 


Boundary  conditions  for  the  bar  free  ends  define  the  wave  number  as 


k]L=nn  (6.8) 

Where  n  =1,2, ...  is  the  number  of  mode. 

The  dispersion  relationship  for  the  flexural  waves  in  thin  elastic  bar  has  the  following 
form  [32]: 


where  E  and  v  are  the  Young’s  and  Poisson’s  moduli  of  the  bar  material  respectively;  p 
is  its  the  density. 


The  Equations  (6.8)  and  (6.9)  allow  to  find  eigenfrequencies  for  flexural  waves 


4L2,j3 


The  longitudinal  strain  du j/dx  produced  by  the  flexural  wave  is  determined  as  [32]: 


d2w 


du  ]/dx  =  z  — —  =  -zk(  A  cos  k{x  cos  corf  . 
dx2 


(6.10) 


Here  z= 0  corresponds  to  the  neutral  plane  of  the  bar. 


We  consider  the  interaction  between  standing  longitudinal  ultrasonic  and  vibration 
flexible  waves.  In  case  of  standing  longitudinal  ultrasonic  wave  its  displacement  is: 


U2~Bcosk2x,  (6.11) 

where  k2=  0)2/ ci  and  c\  is  the  speed  of  the  probing  longitudinal  wave: 

cj-JeTp,  (6.12) 

Substitution  (6.10)  and  (6.1 1)  into  (6)  yields  the  equation  for  the  crack  generated  force,  F, 
at  the  combination  frequencies  (O^  —  +  0)2  : 

F  =  lh  Ac?3  =  alhz0k  ^ k2D2AB cos kxx0  sin  k2x0  cos co3t ,  (6.13) 

We  consider  the  rectangular  crack  passing  the  whole  bar. 


Force  F  applied  to  the  bar  generates  ultrasonic  waves  U3  cosfi)3t  at  the  respective 

combination  frequencies  (Ov  Since  the  wavelengths  of  these  waves  are  greater  than  the 
bar  thickness,  we  assume  that  the  force  F  is  uniformly  distributed  across  the  bar  cross- 
section  at  the  crack  location  x  =  x0  and  generates  plane  (zero-mode)  longitudinal  waves. 


The  resulting  displacement  distribution  along  the  bar  at  the  frequencies  can  be 
presented  in  the  form 


U3= 


Ci  coskpc 
C2  cosk3(L-x) 


if  x<xq 
if  x  >  xq 


(6.1* 


This  form  is  chosen  to  satisfy  the  free  boundary  conditions  at  the  end  of  the  bar 
dUs/dx  |x=o,L  =  0 

and  at  the  cross-section  of  the  applied  force  F  at  x  =  x0.  Coefficients  Cj  and  C2  are 
defined  by  the  conditions  of  displacement  continuity  at  the  cross-section  x  =  x0: 

U3+  =  U3.  (6-15) 

and  the  stress  “jump”  due  to  the  applied  force  F: 

-  03-  =2  F/S  ,  (6.16) 

where  S  is  the  cross-section  area  of  the  tested  plate. 

The  stress-strain  relationship  for  the  longitudinal  deformations  in  the  bar  has  the  form: 

(6.17) 

dx 

The  coefficients  Cj  and  C2  can  be  found.  Substituting  (6.14, 6.17)  into  (6.15, 6.16) 

Finally,  the  amplitude,  U3,  of  the  longitudinal  plane  wave  with  the  combination 
frequencies  C03  is  defined  as: 

AT> 

U3=  z0ak  2k2lhD 2  — — — — —  cos[k3(L  -  x0)]cos(V0)sin(£2x0)  (6-l8) 

1  Sk3  sm{k3L) 

The  term  l/sinfifcjZ,)  in  (6.18)  indicates  resonance  nature  of  the  waves  with  the 
combination  frequencies.  The  dissipation  of  these  waves  can  be  accounted  by  introducing 
complex  wave  number 

k3  =  k  3  (\+i/2Q) ,  (6-19) 

here  Q  is  the  quality  factor  of  the  resonance  oscillations  of  the  bar  at  the  frequencies  (03. 

Besides  this  resonance  term,  the  solution  (6.18)  contains  a  few  more  frequency  dependent 
terms,  which  could  complicate  the  interpretation  of  the  measurements.  This  problem  can 


be  avoided  by  averaging  of  the  measured  nonlinear  response  across  the  frequency  range. 
There  could  be  two  types  of  averaging  yielding  two  criteria  TE  and  TA: 


The  energy  averaging 


(6.20) 


and  the  amplitude  averaging 

N 


SI 


(6.21) 


In  the  experiments  we  used  the  amplitude  averaging  of  the  received  signals  for  the 
following  reason.  Near  the  resonance  the  amplitude  is  proportional  to  the  quality  factor 
Q,  while  the  width  of  the  resonance  curve  is  proportional  to  1  IQ.  Therefore  the  amplitude 
averaging  (integration)  of  the  frequency  response  yields  criteria  independent  (or  weakly 
dependent)  on  the  dissipation  in  the  system. 

As  could  be  seen  from  Eqs.  (6  .18,  6.21)  the  amplitude  of  the  modulation  components  and 
their  frequency  averaged  values  are  linearly  proportional  to  the  amplitude  of  vibration  A, 
amplitude  of  ultrasonic  wave  B,  and  the  crack  depth  /. 

Ta~  U3~  z0ahlD2AB  (6.22) 


These  dependencies  were  verified  U3  (/„ )  experimentally. 

The  side  component  calculation  had  been  conducted  for  several  cases. 

The  polycarbonate  plate 

The  first  example  is  polycarbonate  plate  with  150  mm  length  and  10  mm  thickness.  The 
results  of  side  component  calculation  for  different  crack  positions  are  presented  in  fig.6.2 
-  6.4  (modulation  frequency  1  kHz,  Q=500,  frequency  band  150  -250  kHz)  are 
presented. 

The  variation  of  the  side  components  for  the  frequency  band  300-350  kHz  is  shown  in 
figures  6.5  and  6.6.  It  can  be  noted  when  the  frequency  is  increased  a  side  component 
amplitude  is  decreased  and  the  resonance  peaks  became  smoother 

The  influence  of  tested  object  length  is  illustrated  by  fig.6.7.  This  figure  presents  the 
calculation  results  for  plate  with  250  mm  length.  The  spectra  show  more  resonance  peaks 
in  the  same  frequency  band  then  for  shorter  plate. 
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It  is  observed  that  the  level  of  side  component  depends  on  frequency.  To  eliminate  effect 
of  frequency  on  results  of  measurements  we  used  the  averaging  of  the  side  components 
according  formula  (6.21).  The  averaged  values  of  the  side  components  were  calculated 
for  the  examples  presented  above  and  the  results  are  given  in  figure  6.8.  It  also  can  be 
noted  that  the  averaged  level  is  slightly  depends  on  the  crack  position  and  this  variation 
decreases  with  the  increasing  of  the  frequency  band. 


Figure  6.2.  Effect  of  frequency  on  side  spectral  component  for  different 
crack  positions.  Specimen  length  =150  mm,  frequency  band  from  150  to 
200  kHz,  Q  factor  Q=500. 
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[Figure  6.3.  Effect  of  frequency  on  side  spectral  component  for  different 

crack  positions.  Specimen  length  =150  mm,  frequency  band  from  150  to 
200  kHz,  Q  factor  Q=500. 

■ 
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Figure  6.4.  Effect  of  frequency  on  side  spectral  component  for  different 
crack  positions.  Specimen  length  =150  mm,  frequency  band  from  150  to 
200  kHz,  Q  factor  Q=500. 
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Figure  6.5.  Effect  of  frequency  on  side  spectral  component  for  different 
crack  positions.  Specimen  length  =150  mm,  frequency  band  from  300  to 
350  kHz,  Q  factor  Q=500. 


Figure  6.6.  Effect  of  frequency  on  side  spectral  component  for  different 
crack  positions.  Specimen  length  =150  mm,  frequency  band  from  300  to 
350  kHz,  Q  factor  Q=500. 


Figure  6.7.  Effect  of  frequency  on  side  spectral  component  for  different 
crack  positions.  Specimen  length  =250  mm,  frequency  band  from  150  to 
200  kHz,  Q  factor  Q=500. 


Figure  6.8.  The  dependence  of  average  modulation  level  on  crack  position 
for  two  frequency  ranges.  The  sample  length  is  150  mm. 


6.2.  Acoustic  Crack  Model  for  Scattered  Damage 


The  acoustical  model  of  a  single  crack  was  discussed  in  the  previous  part  of  report.  The 
experimental  results  presented  in  sections  4  proved  that  the  developed  technique  is  also 
applicable  for  detection  of  material  degradation  when  no  cracks  have  been  observed.  The 
model  for  acoustical  properties  of  material  with  a  system  of  microcracks  proposed 
below. 


We  assume  that  tested  specimens  have  uniform  distribution  of  microcracks  (scattered 
damage)  oriented  isotropically  and  orientation  of  microcracks  has  the  same  probability. 


Let’s  first  consider  a  linear  elastic  isotropic  media  (matrix)  with  Young’s  modulus  E0  , 
Poisson’s  ratio  v0  and  the  shear  modulus  G0  =  E0  / 2(1  +  v0).  The  matrix  contains 
cracks,  which,  for  simplicity,  are  assumed  to  be  identical,  circular  of  radius  r  and 
distributed  randomly  (in  the  homogeneous  manner)  relative  to  both  their  locations  and 
orientations.  The  crack  density  is  defined  as 

P  =  ^~  (6-23) 


where  n  is  the  number  of  cracks  in  characteristic  volume  V. 


Based  on  [34]  Young’s  modulus  E  and  Poisson’s  ratio  v 
randomly  oriented  microcracks  with  crack  density  /?is  : 

"'"1  E 


E  =  E, 


16(1-Vq)(1-v0/4) 
9(l  +  vc/2) 


V  =  v„ 


of  the  media  with  system  of 


(6.24) 


The  shear  modulus  G  of  the  cracked  media  is: 
E 


G  = 


(6.25) 


2(1 +  v) 

From  (6.24)  and  (6.25),  it  is  seen  that  the  cracks  transform  the  isotropic  matrix  into  a 
new,  also  isotropic  matrix  with  new  elastic  moduli. 


Our  goal  is  to  find  the  parameter  that  characterizes  interaction  of  ultrasonic  wave 
modulation  with  vibration.  This  nonlinear  acoustic  parameter  y  describes  the  derivation 
“stress-strain”  dependence  of  the  material  from  Hook’s  law.  For  longitudinal 
deformations  this  parameter  is: 

a(s)  =E0  (e+y^/2)  (6.26) 


Therefore,  the  nonlinear  parameter  yean  be  given  as: 
1  d2a 


r  - 


EEs2 


(6.27) 


or  in  terms  of  E  and  a  as: 
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r 


dE_ 

da 


(6.28) 


Under  assumption  that  crack  concentration  and  variation  is  small  the  Eq.(6.24)  can  be 
written  in  form 


E=Eo+AE 

AE=Kp,  (6.29) 

Where 

ri6(l-v02)(l-vn/4)~ 

9(1 +  v,  12) 


The  nonlinear  parameter  can  be  also  expressed  in  term  of  applied  stress  and  crack  density 
on  applied  stress: 

da 

Hence,  if  we  can  find  dependence  of  the  density  of  open  crack  density  on  applied  stress 
a  we  will  be  able  to  determine  nonlinear  parameter  y. 

Let’s  consider  bilinear  model  approach.  According  this  model  all  cracks  are  opened 
under  the  tension  stress  (a> 0)  and  are  closed  under  the  compression  stress  (cr<0).  This 
model  cannot  adequately  describe  observed  phenomena  and  we  will  consider  more 
realistic  model  that  includes  the  presence  of  residual  stresses.  Suppose  that  there  are 
residual  stresses  in  material  and  crack  can  be  opened  if  the  applied  stress  exceeds  residual 
stress.  Let  describe  residual  stress  by  function  of  stress  distribution  W(a*). 

It  means  that  W(a*)da*  is  a  probability  that  residual  stress  in  any  crack  is  in  the  interval 
from  a » to  a*+  da*. 

Let  consider  the  Normal  distribution  of  residual  stresses: 


W(a.)  = 


la: 


2\ 


(6.31) 


Where  A  is  variance  of  residual  stress  distribution.  Suppose  that  the  average  value  of  the 
residual  stress  is  close  to  zero. 

First  let’s  consider  the  case  when  the  longitudinal  stress  a  is  applied  to  the  thick  bar 
having  uniform  crack  distribution.  Let  consider  the  cracks  having  the  normal  n  oriented 


at  an  angle  cp  to  the  z  axis  (Fig.6.9).  The  normal  stress  applied  to  these  cracks  is  defined 
by  equation: 

<rm  -  a  cos2  <P 


Figure  6.9.  The  coordinate  system  described  crack  orientation 


The  crack  density  having  orientation  of  the  normal  between  (cp,0)  and  (cp+dcp,  0+d0  )  is 
defined  by  equation: 

N(<p,e)=  (6.32) 

An 

(0<cp<7t,  O<0<27l). 

The  crack  is  opened  if  the  normal  stress  applied  to  the  crack  exceeds  residual  stress  acted 
on  the  crack.  The  density  of  the  opened  cracks  can  calculated  as 

2  it  2x 

p(cr)  =  Po  £C°S  *  W(a.  )dcr ,  J  \H((p,6)dddL<p  (6.33) 

0  0 

Where  p0  is  the  crack  density  when  all  cracks  are  opened. 

The  nonlinear  parameter  y  can  be  found  by  derivative  of  this  equation  according  (6.30) 
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=  k^B-- 


r  =  K 


p°  r r w  ^■'>da-  jsin«*^ 


da 


Ida 


(6.32) 


This  differentiation  with  respect  to  a  parameter  can  be  conducted  by  Leibnitz’s  rule 


a  v(A)  »(A) 

—  |/(X,2)^=  J 


df(x,X) 

dl 


dx  +  f(v,X)— 

oA 


(6.33) 


Rearranging  (6.33) : 

It 

y  =  —  jKTF (0)  Jcos 2  p  sin<pd<p  =  —  ATT (0)  (6.34) 

Substitution  of  residual  stress  distribution  from  (6.31)  gives  the  final  answer  for  nonlinear 
parameter  of  the  cracked  media: 


(6.35) 


It  can  be  noted  that  the  nonlinear  parameter  is  proportional  to  the  crack  density  and 
inverse  proportional  to  the  variance  of  residual  stress  distribution. 


7.  Impulse  Nonlinear  Method  of  Crack  Location 


The  NEWS  technique  used  CW  ultrasonic  signal  and  described  above  proved  high 
sensitivity  for  “pass  -fail”  test  but  cannot  be  used  to  locate  a  crack  position.  We  propose 
that  the  crack  location  can  be  defined  by  a  pulse  modulation  method.  This  technique 
provides  precise  defect  location  and  the  concept  is  shown  in  fig.7.1. 


Vibrator 


Crack 


Signal  from 
hole 


Signal  from 
crack 


Received  signal 


t 


Figure  7.1 .  Schematic  of  nonlinear  acoustic  technique  for  flow  location. 

The  acoustical  transducer  in  conventional  technique  produces  an  ultrasound  impulse 
propagating  in  the  specimen.  This  impulse  is  reflected  from  different  kinds  of 
inhomogeneities.  For  example,  and  the  reflected  signals  from  a  hole  and  a  crack  are 
presented  in  fig.7.1.  Reflected  signals  came  at  different  time  but  there  is  no  way  to 
distinguish  what  signal  is  reflected  from  the  crack  and  what  signal  from  the  hole.  These 
signals  can  be  easily  distinguished  if  mechanical  vibration  can  be  applied  to  the 
specimen.  Due  to  vibration  the  amplitude  and  phase  of  the  ultrasound  signal  reflected 
from  crack  will  be  modulated,  while  the  signal  reflected  from  linear  defect  (the  hole)  has 
no  modulation.  This  situation  is  illustrated  in  fig.7.2. 
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Figure  7.2..  The  time  tracks  of  the  ultrasound  signal  reflected  at  different  phases  of  the 
vibration. 

The  upper  part  of  this  figure  presents  the  harmonic  vibration  with  low  frequency.  The 
lower  part  of  the  fig.  7.2  shows  signals  reflected  from  the  defects  at  different  times.  The 
first  pair  is  reflected  at  the  moment  when  the  applied  vibration  stress  reaches  the 
maximum.  At  this  time  the  crack  is  compressed  and  the  signal  reflected  from  the  crack 
has  minimal  amplitude.  The  second  pair  of  impulses  is  shown  at  time  when  the  vibration 
stress  is  small  and  there  is  no  change  in  the  reflected  impulse.  The  third  pair  is  extension 
phase  of  the  vibration,  at  that  moment  the  crack  is  maximally  opened  and  the  reflection 
signal  is  higher  than  in  previous  phases.  The  modulation  of  this  signal  manifests  the  crack 
presence  allows  us  to  distinguish  the  crack-like  defects  from  other  inhomogeneities. 

The  test  demonstrated  this  approach  was  carried  out  by  using  a  piezoceramic  transmitter 
of  ultrasound  wave. 

First  experiments  were  conducted  on  cylindrical  aluminum  rods  of  (length  =  91cm; 
diameter  1  cm).  The  rods  with  a  crack  and  with  a  notch  were  investigated. 

The  experimental  setup  is  shown  in  fig.  7.3.  The  high  frequency  piezoceramic  transducer 
and  receiver  were  used  for  radiation  and  reception  of  the  ultrasound  wave  with  frequency 
of  1.2  MHz.  Low-frequency  vibrations  of  the  rod  were  excited  by  a  vibrator  ( mass  557 
g).  The  vibrator  can  provide  different  excitation  frequencies  and  for  our  experiment  the 
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frequency  around  10.5  Hz  was  chosen,  which  is  the  main  flexural  resonance  frequency 
of  the  rod  loaded  with  the  mass  of  vibrator  (557  g). 


Figure  7.3.  The  experimental  setup  to  impulse  nonlinear  method  of  crack  location. 


The  time  tracks  of  the  received  signals  reflected  from  the  notch  and  the  crack  are 
presented  in  fig.7.4.  It  is  seen  that  there  is  no  clear  difference  between  the  signals 
reflected  from  the  defects. 
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Figure  7.4.  The  time  tracks  of  the  signals  reflected  from  the  notch  (a),  and  the 
crack  (b) 


The  received  ultrasound  signals  passed  trough  the  gating  device,  which  allows  the  receipt 
of  signals  during  a  definite  time  interval.  In  the  experiment  this  time  interval  was  chosen 
as  6  ms.  The  starting  time  can  vary  to  allow  us  to  pick  up  the  signal  from  different 
distances.  From  the  gating  device  the  signal  was  sent  to  a  peak  detector  and  the  sequence 
of  the  impulses  was  observed  after  detection.  This  sequence  is  modulated  by  low 
frequency  vibration  if  a  signal  is  coming  from  the  crack.  The  manifestation  of  such 
modulation  can  be  easily  seen  in  the  spectrum  of  the  detected  signal.  The  presence  of 
modulation  gives  a  harmonic  component  with  the  vibration  frequency  while  there  are  no 
such  components  in  the  spectrum  of  the  signal  reflected  from  the  hole  ( see  Figure  7.5.). 

It  is  clear  that  this  effect  can  be  used  to  distinguish  cracks  from  other  kinds  of 
inhomogeneities. 


frequency,  Hz 


Figure  7.5a.  The  spectrum  of  the  signal  reflected  from  the  crack. 


frequency,  Hz 


Figure  7.5.b.  The  spectrum  of  the  signal  reflected  from  the  hole. 


The  another  experiments  have  been  conducted  by  using  higher  probe  ultrasonic 
frequency  (3  MHz)  that  allowed  to  increase  the  spatial  impulse  resolution.  The  steel  plate 
having  length  about  145  mm  and  width  60  mm  was  used  for  demonstration.  The  crack  in 
the  plate  was  grown  in  UIC  under  cyclic  loading.  The  sensor  was  placed  on  the  side  of 
plate  and  the  technique  described  above  was  applied  to  get  received  impulse. 

Figures  7.6  and  7.7.  demonstrates  the  result  of  measurements.  The  linear  reflection 
signals  and  side  level  in  nonlinear  processing  are  presented  in  3D  form.  Axis  Z  is 
transverse  coordinate  of  plate  where  sensor  was  placed.  The  distance  between  sensor  and 
point  of  reflection  is  measured  along  X  axis.  This  coordinate  is  connected  with  time  of 
the  reflected  signal  by  using  formula: 
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x=2ct 

Where  c  is  the  sound  speed  of  the  tested  signal.. 

It  is  clear  seen  that  the  linear  technique  gives  very  similar  signal  s  from  the  crack  and  the 
hole  whereas  the  nonlinear  processing  provides  very  distinguish  difference  between 
cracks  and  defect  without  nonlinear  properties  (like  holes). 

The  amplitude  dependence  for  nonlinear  signal  was  also  checked.  The  fig.7.8.  presents 
the  result  nonlinear  imaging  for  two  different  level  of  low  frequency  vibration.  It  is  seen 
that  the  side  component  value  increases  with  vibration  amplitude. 

The  fig.7.9.  gives  the  dependence  of  the  side  component  signal  versus  vibration 
displacement.  It  is  seen  that  this  dependence  is  close  to  the  linear  one.  It  confirms  the 
theoretical  model  presented  in  section  6  of  the  report . 
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Figure  7.6.  Comparison  of  linear  and  nonlinear  location.  Hole  is  situated  closer  to  the 
sensor  then  the  crack.  The  nonlinear  signal  from  the  crack  significantly  exceeds  the 
nonlinear  signal  from  the  hole.  The  signals  near  the  plate  end  are  due  to  the  reflection 
from  the  end. 


Figure  7.7.  Comparison  of  linear  and  nonlinear  location.  The  crack  is  situated  closer  to 
the  sensor  then  the  hole.  The  nonlinear  signal  from  the  crack  significantly  exceeds  the 


nonlinear  signal  from  the  hole.  The  signals  near  the  plate  end  are  due  to  the  reflection 


from  die  end. 


Figure  7.9.  The  dependence  of  the  side  band  signal  from  the  vibration  amplitude. 
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8.  Kinematics  and  Mechanical  Design  of  the  Instrumental 
Hammer 

This  part  of  report  gives  the  calculation  of  the  instrumental  hammer  parameter  that  can 
provide  necessary  impact  excitation 

The  hammer  has  to  provide  the  following  parameters: 

Impact  Impulse:  from  200  to  2000  N  msec 

Distance  to  sample:  from  5  to  20mm 

Source  of  energy:  spring. 

8.2.  Process  of  impact 

The  main  parameter  that  the  hammer  has  to  provide  is  the  impact  impulse 
S  =  F  t. 

Where 

F  is  force,  t  is  the  impulse  duration. 

Because  the  hummer  head  is  bounced  after  the  impact 

S  =  m  (u+v), 
where 

m  -  mass  of  the  hammer  head 
u  -  speed  after  impact 
v  -  speed  before  impact 

Coefficient  of  restitution : 


k  =  u/v 

When  the  speed  is  v  =3  m/s ,  k  =  0.6-0.8  for  steel  -  steel  impact 


u  =kv 


S  =  m  (k+l)  v 

Speed  of  the  hammer  head  to  provide  required  impulse  is 
v=S  /m  (k+l)  (8.1) 


8.3.  Calculation  of  spring  parameters. 


Lo  =  free  length  of  spring 

Xo  =  deformation  of  spring  before  impulse 

Change  of  kinetic  energy  of  system 

mv"/2  -  mv"o/2  =  SA  ,  where 

vo  =  initial  speed 
v  =  final  speed 
A  =  total  work 

Work  of  spring  force  F: 


F  =  CX, 


C  =  stiffness 
X  =  deformation 

A  =  C  (X"o-X"l)  /  2 

At  the  end  of  release  deformation  of  spring  is  zero  and  the  work 
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A  =  C  X"o  /  2 


Therefore, 

mv"  /  2  =  C  X"o  /  2 

Speed  at  the  end  of  flight: 

v  =  Xo  Vc/in 

From  (  8.1)  and  (8.2) 

S  /  m  (K+l)  =  Xo  VC/m,  where 
Stiffness: 

C  =  [S  /  Xo  (K+l)]"  /  m 
If: 

S  =  0.5N  Sec 

Xo  =  0.04  m  (1.6”) 

K=0.7 
m  =  0.1  kg 
C  =  550  N/m 

From  (8.2): 

v  =  2.83  m/sec 


Work  is  made  by  spring  force  and  weight  of  the  worker  (shock 
from  up  to  down) 

SA  =  C  X"o  /  2  +  m  g  Xo  =  mv"  /  2, 


where 


v  =  V(CX"o/m+2gXo) 

If: 

m  =  0.1  kg 
C  =  500  N/m 
Xo  =  0.04  m 
V  =  2.96  m/sec 

Difference  is  ( 2.96  -  2.83  )  /  2.83  =  0.047  or  4.7% 


85 


Conclusion:  weight  doesn’t  significantly  affect  the  speed 

8.4.  Shock  worker  without  free  flight 

Range  of  regulation  by  one  spring. 

Schematic  is  shown  in  Figure  8.2. 

Lo  =  free  length  of  spring 
AX  =  motion  of  shock  worker 

X  =  deformation  of  spring  for  regulation 

Xo  =  deformation  of  spring  before  impact 

XI  =  deformation  of  spring  after  impact 

51  =  min.  impulse  from  the  spring 

52  =  max.  impulse  from  the  spring 


Figure  8.2.  Regulation  by  one  spring 

For  SI:  Xo  =  AX 
XI  =0 

For  S2:  Xo  =  X  +  AX 
XI  =X 
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Change  ofkinematical  energy  of  system 
mv72  =  C  (X"o-X"l)  /  2,  where 

v  =  VC(X"o-X"l)/m  (8.4) 

V2/Y1  =  S2/  SI  =Ks,  where 
VI  and  V2-min.  and  max.  speed 

V2  =  Ks  VI  (8.5) 

For  min  SI: 

VI  =  <  C(AX"-0)/m 
For  max  S2 

V2  =  V  C(X+AX)"-X")/m  =  VC(AX"+2AX  X)/m 
V2  =  Ks  VI 

VC(AX"+2AX  X)/m=  Ks  AX,  where 
X  =  deformation  of  spring  for  regulation 
X  =  (Ks"-1)  AX/2 
For  ratio: 

Ks  =  S2/Sl  =2000/200=  10 
And  AX  =  25mm 

X  =  1 1 12mm.  That  is  too  much 

For  ratio: 

Ks  =  1.5 

And  AX  =  25mm 
X  =  15mm 
For  ratio: 

Ks  =2  X=37.5mm 

Therefore  one  spring  can  regulate  for  Ks  in  the  range  from  1.5  to  2 


Variation  of  speed  by  change  of  the  distance  between  worker  and 
sample  (internal  range  of  spring  deformation) 


Figure  8.3.  Change  of  the  distance  between  shock  worker  and 
sample  (internal  range  of  spring  deformation) 

Lo  -  free  length  of  spring 
AX  -motion  of  shock  worker 

X  -  deformation  of  spring  for  regulation 
Ls  -distance  to  sample 

AXs  -unused  deformation  of  spring 

Mark  in: 
i  =  Ls  /  AX 

Ls  =  AXi 
AXs  =  AX  -  AX  i 

From  4: 

V  =  V  C(Xo"-Xl")/m 

Xo  =  deformation  of  spring  before  impact,  Xo  =  AX  +  X 

XI  =  deformation  of  spring  after  impact,  XI  =  AXs+X  =  AX-AX  i  +X 


Markin  A  =  AX+X, then 


Below  is  the  speed  of  worker  at  the  instant  of  impact  when 
distance  between  worker  and  sample  is  changed. 


=  500  N/m, 

3 

ii 

p 

P 

AX  = 

0.03  m, 

X  =  0.01  m, 

A  =0.04m 

i 

0 

0.17 

0.4 

0.67 

0.8 

Ls 

0 

5 

12 

20 

24 

V 

0 

1.37 

2.02 

2.45 

2.59  2 

Note  that  speed  varies  significantly  when  the  distance  between 
worker  and  sample  is  changed  (for  the  distance  from  5  to  20  mm 
the  difference  in  speed  is  44%). 
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In  free  flight  the  speed  of  the  worker  varies  insignificantly  when  the 
distance  between  the  worker  and  sample  is  changed. 

Friction  force: 

F  =  m  g  f,  where: 

m  =  mass  of  shock  worker 
g  =  acceleration  of  gravity 
f  =  coefficient  of  friction 

F=m  a,  where  a  =  acceleration  and  a  =  g  f 

Speed  of  flight  is  defined  by : 

V"  =  2a  L,  where  L  =  distance 

V"  =  2g  f  L 

After  differentiation  we  have: 

2  V  AV=2gf  AL 
and  speed  change  is 
A V  =  g  f  AL  /  V 

If:  g  =  9.8  m/sec",  f  =  0.1,  V  =  2.8  m/sec,  AL  =  0.04  m 

AV  =  0.014  m/sec  and 
AV/V  =  0.005  =  0.5% 

In  free  flight  speed  and  impulse  are  changed  insignificantly 
Shock  worker  time  of  flight  is: 
t  =2  L  /  V 

If:  L  =  0.04m,  V  =  2.8m/sec,  t  =  0.028sec 
Definitions  of  common  parameters: 

Min.  impulse:  min  SI 
Deformation  of  spring  before  shock 


Xo=Xl 
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Max.  impulse  max  S2 
Deformation  of  spring  before  shock 
Xo  =  X  2  =  X  +  XI 
Formula  (8.2)  gives: 

VI  =  XTTC/m 
V2  =  (X  +  XTjVC/m 
It  is  followed  from  formula  (8.5): 

V2  =  Ks  VI 

V2  =  (X  +  XI)  a/C  /  m  =Ks  XlVC/m ,  where 
X  +  XI  =  Ks  XI 

Deformation  of  spring  to  regulate  impulse  from  min  SI  to  max  S2  : 

X  =  XI  (Ks-1)  (8.8) 

If:  XI  =  15  mm 
Ks  =1.8, 

X  =  12mm 
X2  =  27mm 


Ranges  of  regulations  (Ks  =  1.8) 
4  springs  are  used  for  regulation 


#  spring  min  SI  max  S2 

1  0.2  0.36 

2  0.36  0.65 

3  0.65  1.17 

4  1.11  2 

Stiffness  is  defined  by  formula  (3): 


C  =  [S/Xo(K+l)]7m 
Mass  of  worker:  m  =  0.1kg 


Force  of  spring: 

PI  =CX1;  P2  =  CX2 


Values  of  common  parameters  : 

#  Spring  SI, 2  N  sec  XI, 2  mm  C  N/m  PI, 2  kg 


0.2 

15 

615 

0.5 

0.36 

27 

615 

1.7 

0.36 

15 

2000 

3 

0.65 

27 

2000 

5.4 

0.65 

15 

6500 

9.8 

1.17 

27 

6500 

17.5 

1.11 

15 

19000 

28.4 

2 

27 

19000 

51.1 

From  the  drawing  of  monoblock  in  figure  8.5: 

Spring  directly  reacts  on  the  worker.  Lift  of  worker  can  be  done  manually  or  by 
electromagnet.  Worker  is  fixed  in  top  by  lock.  When  electromagnet  goes  down  it  releases 
lock  and  shock  worker  moves  down  to  sample  fast.  The  lock  can  be  released  manually. 
Electromagnet  is  turning  on  at  once  right  after  impact  and  lifts  shock  worker  up. 

8.6.  Modification  with  free  flight  of  worker  with  cock. 

Device  has  2  elements  on  one  axis:  worker  and  cock.  Compressing  spring  directly  reacts 
on  the  cock.  Lift  of  cock  is  done  manually  or  by  electromagnet.  The  cock  is  fixed  by 
lock  on  top. 

When  electromagnet  goes  down  it  releases  lock  and  cock  moves  down  fast  to  shock 
worker .  The  lock  can  be  released  manually. 

Electromagnet  is  turning  on  at  once  after  shock  and  lifts  cock  up.  Shock  worker  goes  up 
after  impact  and  is  held  up  by  light  spring.  To  avoid  going  down,  shock  worker  has  soft 
pad.  This  device  has  less  axial  dimension. 

Design  is  given  in  Figure  8.6. 
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8.7.  Modification  with  free  flight  of  worker  with  rotating  arm. 

Device  has  2  elements:  worker  and  rotating  arm .  Extension  spring  directly  reacts  on  the 
arm.  Lift  of  arm  is  done  manually  or  by  electromagnet.  The  arm  is  fixed  by  lock  on  top. 
When  electromagnet  goes  down  it  releases  lock  and  arm  moves  down  fast  to  shock 
worker .  The  lock  can  be  released  manually. 

Electromagnet  is  turning  on  at  once  after  shock  and  lifts  arm  up. 

Shock  worker  goes  up  after  impact  and  is  held  up  by  light  spring.  To  avoid  going  down, 
shock  worker  has  soft  pad.  This  device  has  less  axial  dimension. 

Design  is  given  in  Figure  8.7. 


Figure  8.5.  Monoblock 


Figure  8.7.  Modification:  free  flight  of  shock  worker  with  rotating 
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9.  Conclusion. 

The  successful  fulfillment  of  Phase  I  project  allows  to  propose  Phase  II  that  will  provide 
an  advanced  development  of  new  technology  and  instrumentation  for  reliability 
assessment  As  demonstrated  in  Phase  I  the  concepts  are  feasible  but  need  significant 
refinement 'and  optimization  to  bring  them  to  a  commercially  usable  state.  When  brought 
to  this  state  at  the  end  of  Phase  II  the  suggested  technique  will  find  wide  use  for  life  time 
prediction  and  reliability  assessment.  It  is  our  belief  that  nonlinear  methods  may  be 
routinely  used  in  applications  as  diverse  as  monitoring  reactor  containment  wall  damage; 
inspecting  aircraft  and  spacecraft;  observing  fatigue  damage  in  buildings,  bndges, 
tunnels,  gas  and  oil  pipe  lines;  et  cetera. 

The  information  formulated  during  Phase  II  will  allow  transfer  of  technology  to  other  Air 
Force  and  military  entities  and  to  the  commercial  market.  In  Phase  III,  Smart  Material 
Design  plans  to  market  different  kinds  of  NEWS  devices  started  from  small  and  lower 
power  devices  for  single  part  testing  up  to  large  and  powerful  equipment  infrastructural 
testing  (including  large  parts  of  aircrafts,  bridges,  buildings,  etc).  Phase  II  proposal  will 
establish  a  solid  foundation  for  successful  commercialization. 

There  are  potentially  a  large  number  of  applications  of  enormous  economic  and  safety 
impact  that  will  evolve  from  non-linear  acoustic  research  proposed  m  this  project.  It  is 
not  understatement  to  suggest  that  within  10  years  non-linear  methods  for  the  reliability 
assessment  of  aging  structures  may  be  routinely  used  in  applications  as  diverse  as 
inspecting  aircraft  and  spacecraft,  monitoring  reactor  containment  wall  damage, 
observing  fatigue  damage  in  buildings,  bridges,  tunnels,  gas  and  oil  pipe  lines,  etc. 

The  proposed  nonlinear  acoustic  approach  is  much  more  sensitive  to  fatigue  damage  than 
the  commonly  used  ultrasonic  non-destructive  techniques.  The  sensitivity  of  non-linear 
wave  methods  to  the  appearance  and  progression  of  material  damage  is  orders  of 
magnitude  higher  than  that  of  conventional  methods  of  non-destructive  testing.  ^Thus, 
measurement  of  this  non-linear  behavior  may  well  be  the  most  sensitive  method  for  toe 
study  of  fatigue  damage  initiation  and  cracking.  That  is  why  we  are  confident  that 
national  research  laboratories  and  universities  will  be  also  interested  in  implementation  of 
proposed  non-linear  acoustic  technique  for  the  fundamental  research  projects. 

Preliminary  investigations  with  non-linear  acoustic  approach  for  toe  damage  assessment 
were  conducted  on  request  and  with  cooperation  of:  Boeing;  Shell;  British  Airspace, 
Dynamic  Response  Systems,  Incorporated;  Ford;  Los  Alamos  National  Laboratory;  Gas 
Research  Institute;  Purdue  University  Calumet;  Electric  Power  Research  Institute; 
University  of  Illinois  at  Chicago;  and  toe  National  Institute  of  Standards  and  Technology. 
Contacts  have  been  initiated  with  Wright-Patterson  Air  Force  Base;  American  Petroleum 
Institute;  LaserSonix  Technologies,  Incorporated;  Iowa  State  University,  and  Texaco. 
The  positive  results  of  preliminary  investigations  are  well  above  expectations  and  all 
above  cited  companies  support  the  project  and  expressed  commitments  to  toe  use  of  non¬ 
linear  acoustic  instrumentation  when  it  will  be  available  on  toe  market.  Dynamic 
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Response  Systems,  Inc.  and  LaserSonix  Technologies,  Incorporated,  expressed  then- 
commitments  to  participate  in  manufacturing  and  marketing  of  the  instrumentation. 
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Appendix 


Support  Letters 


Wayne  E.  Woodmans*! 
Technical  Fellow 
Operations  Technology 


Boeing  Commercial  Airplane  Group 
P.O.  Box  3707  MS  73-30 
Seattle,  WA  96124-2207 


8  February  1099 
Dr.  Yuiian  Kin 

Professor  of  Mechanical  engineering 
Purdue  Calumet  University 
Hammond,  IN  46323 

Dear  Professor  Kin, 

Professor  Alexander  Sutin  has  informed  me  that  he  will  be  working  with 
your  on  a  proposal  to  foe  NSF  dealing  with  nonlinear  Acoustic  NDE.  I  have 
supplied  adhesively  bonded  test  articles  to  Prof.  Sutin  for  evaluation  wWi  his 
non-linear  acoustic  techniques,  and  I  am  very  interested  in  the  results  he 

has  obtained,  (believe  foe  proposed  work  will  use  nonlinear  acourac 

techniques  for  foe  detection  and  investigation  of  fatigue  damage.  This  is  an 
area  in  which  Boeing  is  very  interested,  particularly  for  the  in-service 
detection  of  fatigue  damage  to  aircraft  components.  Measurements^ 
cumulative  fatigue  damage  or  predictions  of  the  remaining  lifetime  of  a 
component  or  structure  would  be  very  valuable  to  the  aircraft  industry. 

I  was  impressed  with  foe  results  Professor  Sutin  has  obtained  with  Ms 
nonlinear  techniques,  and  I  look  forward  to  foe  outcome  of  your  application 
of  nonlinear  techniques  to  investigations  of  fatigue  damage. 


Sincerely  yours, 

WayneWoodmansee 
Sr.  Technical  Fellow 
426-965-3188 

wayne.e.woodmansee@boeing  .com 


Honeywell 


Engines  &  Systems 
Honeywell 

111S.  34th  St,  P.O.  Box  521  SI 
Phoenix,  A Z  85072-2181 


Dr.  Yulian  B.Kin 
President 

Smart  Material  Design,  Inc. 

31 14  Belwood  Lane, 

Glenview,  Illinois  60025 

fax:  847  729  1813 

e-mail:  Vin@calumet.purdue.edu 


April  18,2000 

RE:  Quantitative  Non-Destructive  Evaluation  and  Reliability  Assessment  of  the 
Aging  Aircraft  Structure  Components 
Principal  Investigator:  Dr.  Alexander  Sutin 

Dear  Dr.  Kin, 

I  am  writing  this  letter  in  support  of  the  further  development  of  above  citied  project.  The 
suggested  technique  appears  to  have  significant  potential  for  nondestructive  inspection  of 
aircraft  engine  parts.  If  proven  successful,  we  expect  that  this  technique  may  open  a  new 
opportunity  in  fatigue  crack  detection  applicable  to  Honeywell  components. 

Where  appropriate,  Honeywell  will  provide  aircraft  engine  parts  and  crack  samples  (on 
loan)  for  testing/demonstration,  and  where  possible,  we  will  also  provide  engineering 
support  of  the  Nonlinear  Acoustic  technique  in  comparison  assessments  with  existing 
techniques  and  demonstration  on  applicable  engine  hardware. 


Sincerelv_yours 

42a 


Tim  Duffy  1|\ 

Manager,  Advanced  Materials  Technologies 
Materials  &Process  Engineering 
Honeywell  Engines  and  Systems 
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Nonlinear  Elasticity  Project 

Group  EES-4 

Earth  and  Environmental 

Sciences  Division 

Los  Alamos  National  Laboratory 

Los  Alamos,  New  Mexico  87545 
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Doctor  Yulian  B.  Kin 
Purdue  University  Calumet 
Engineering  Department 
Hammond,  IN  46323 
Phone:  219  989  2684 
e-mail:  inn@cai™mftt  purdue.edu 

RE:  DEVELOPMENT  OF  NONLINEAR  ACOUSTIC  INSTRUMENTATION 


Dear  Dr.  Kin, 

I  am  writing  in  support  of  the  further  development  of  the  above-cited  project.  I  understand 
that  Smart  Material  Design,  Inc.,  Purdue  University,  and  University  of  Illinois  are 
intending  to  extend  their  research  in  nondestructive  characterization  of  materials  and  aging 
structures.  Los  Alamos  National  Laboratory  is  one  of  the  few  laboratories  in  the  world 
working  in  this  direction  as  welL  We  see  countless  applications  to  aging  structures, 
quality  control  in  assembly  lines,  damage  inspection  from  earthquakes,  etc.,  and  see 
collaboration  as  a  benefit  to  afl. 

We  have  had  a  long-term  collaboration  with  Dr.Sutin.  He  has  conducted  experimental 
research  on  new  methods  for  nonlinear  acoustic  Nondestructive  Testing  in  Los  Alamos 
and  elsewhere,  and  has  been  the  primary  person  in  the  world  responsible  for  development 
of  wave  modulation  methods.  As  you  know,  we  have  also  been  working  with  Dynamic 


Resonance  Systems,  Inc.,  with  the  intention  of  ultimately  marketing  a  device.  May  I 
suggest  that  Los  Alamos,  DRS,  U.  Illinois  Chicago,  Purdue  and  your  group  cooperate  on 
this  venture? 

I  would  be  happy  to  discuss  details  with  you  in  person  or  by  telephone. 


Sincerely, 


Paul  Johnson 
Leader 

Nonlinear  Elasticity  Project 
paj@lanl.gov 


Equipment  available  at  LANL 

2  DRS  Resonance  Ultrasound  Spectroscopy  units 

4  ENI  high  power,  high  frequency  amplifiers 

2  Crown  audio  amplifiers 

Hafler  audio  amplifier 

Krohn-Hite  audio  and  ultrasonic  power  amp 

5-10  wave  source  generators  (HP,  Analogic,  Krohn-Hite) 

Analogic  16  bit  digital  oscilloscopes 

2  LeCroy  8  bit  digital  oscilloscopes 
several  Tektronix  analog  scopes 
Tectronix  waveform  analyzer 

HP  Impedence  Analyzer 

Analogic  arbitrary  function  generator 

3  lock-in  amplifiers 
HP  pulse  generator 

Bruel  and  Kjaer  Accelerometers,  pre-amps  and  amplifiers. 

Three  Macintosh  G3’s  for  experimental  control  with  Lab  View  control  software. 
2  complete  Nonlinear  Resonant  Ultrasound  Spectroscopy  (NRUS)  experimental 

ups 

1  Nonlinear  Wave  Modulation  experimental  set-up. 

Numerous  piezoelectric  ceramics,  connectors,  etc.  for  experiments. 

Numerous  samples  of  damaged  and  undamaged  materials. 

Vacuum  system  and  supplies  for  experiment. 

Numerous  passive  and  active  wave  filters 
signal  processing  software  (MatLab  and  IGOR) 

Access  to  an  Origin  2000  (100  parallel  processing  machine)  for  modeling 


united  states  department  of  commerce 

National  Institute  of  Standards  and  Technology 

325  Broadway 

Boulder,  Colorado  80303-332B 


Reply  to  the  attention  of: 

January  21, 1999 
Dr.  Yulian  Kin 

Professor  of  Mechanical  Engineering 
Purdue  Calumet  University 
Hammond,  IN  46323 


Dear  Prof.  Kin: 

r  am  writing  with  regard  to  Dr.  Alexander  Sutin  of  the  Stevens  Institute  of  Technology 
LThTp  Jcipaion  to  the  NSF  proposal  entitled  “Nonltoear  Aconstic  Instrumentation 
for  the  Detection  and  Investigation  of  Fatigue  Damage.” 

In  mv  position  as  a  physicist  at  the  National  Institute  of  Standards  &  Technology,  l  am 
developing  new  measurement  techniques  in  the  field  of  nonlinear  acoustics.  My  researc 
is  aimed  primarily  at  using  these  techniques  for  the  nondestructive  charactenzanon  of 
materials  Dr.  Sutin  has  made  several  notable  contributions  to  this  same  field.  For 
instance,  he  has  demonstrated  the  use  of  the  vibroacoustic  method  for  detecting  feugue 

damage. 

Although  I  have  been  aware  of  his  work  since  I  began  working  in  nonlinear  acoustics  in 
1995  I  did  not  meet  Dr.  Sutin  until  1997.  Since  then,  we  have  maintained  contact 
through  personal  correspondence  and  meetings  at  professional  conferences.  Based  on  my 
interactions  with  Dr.  Sutin,  he  seems  to  be  an  excellent  researcher,  with  strong  skills  m 
hnth  the  experimental  and  theoretical  areas.  He  appears  eager  to  meet  and  overcome 
technical  challenges,  and  is  very  enthusiastic  about  realizing  the  full  potential  of 
£ 52222  for  materials  characterization.  Therefore,  in  this  letter  I  wishto  voice 
mv  ciinnort  for  his  part  in  the  proposal.  I  am  confident  that  his  performance  will  be 
excellent,  and  that  he  will  work  towards  accomplishing  the  objectives  to  the  very  best  o 

his  ability. 


Sincerely  yours, 

Donna  C.  Hurley,  Ph.D. 

Materials  Reliability  Division 
National  Institute  of  Standards  &  Technology 


I  FORD 

George  Mozurkewich 
Ford  Research  Laboratory 
P.  O.  Box  2053,  Mail  Drop  3028 
Dearborn,  MI  48121-2053 

January  22, 1999 

Dr.  Yulian  Kin 

Professor  of  Mechanical  Engineering 
Purdue  University  Calumet 
2200  169*  Street 
Hammond,  IN  46323-2094 

RE-  NSF  Proposal  entitled  “NONLINEAR  ACOUSTIC  INSTRUMENTATION 
FOR  THE  DETECTION  AND  INVESTIGATION  OF  FATIGUE  DAMAGE” 

I  Dear  Dr.  Kin, 

I  am  writing  in  support  of  the  proposal  cited  above.  Because  of  my  familiarity  with  Dr. 
Sutin's  work,  I  arranged  for  Ford  Motor  Company  to  provide  him  and  his  colleagues  with 
|  engine  heads,  composite  plates,  and  other  parts  for  preliminary  experiments.  The  results 
of  those  tests  were  very  impressive. 

I  Consequently,  I  anticipate  that  this  work  will  have  significant  scientific  impact  in  the 
areas  of  metallurgy,  materials  science,  and  nondestructive  evaluation  (NDE). 
Furthermore,  I  expect  that  further  development  of  these  techniques  will  in  due  course 
j  become  of  great  economic  importance  to  industry. 

As  a  scientist  in  the  Ford  Research  Laboratory,  I  plan  to  maintain  close  contact  with  Dr. 

I  Sutlin  during  the  progression  and  completion  of  this  work,  and  I  hope  that  his  techniques 
■  will  also  prove  useful  in  my  own  endeavors. 

|  Sincerely, 

I 


Intemet:gmozurke@ford.com 

Tel:  (313)  845-5038 

Fax  (313-594-6863 


George  Mozurkewich 
Senior  Technical  Specialist 
(313)845-5038 


Office  of  the  Chancellor 


December  9,  1999 


Dr.  Spencer  T.  Wu 
Program  Manager 

Air  Force  Office  of  Scientific  Research 
801  N.  Randolph  Street,  Room  732 
Arlington  VA  22203-1977 

Dear  Dr.  Wu: 

At  the  request  of  Dr.  Yulian  Kin,  Professor  of  Mechanical  Engineering  at  Purdue 
University  Calumet,  I  am  informing  you  of  the  University's  support  of  the  following 
project:  Quantitative  Non-Destructive  Evaluation  and  Reliability  Assessment  of  Aging 
Aircraft  Structure  Components.  Professor  Kin  is  a  respected  member  of  our  faculty 
and  maintains  a  high  level  of  research  activity  that  has  been  beneficial  to  him,  to  the 
university,  and  particularly  to  our  students. 

My  review  of  the  proposal  gives  me  confidence  that  control  instrumentation  based  on  a 
non-linear  acoustical  approach  will  be  widely  used  for  quality  control  in  basic 
university  research,  in  industrial  applications,  and  in  fatigue  damage  assessment  in 
aoing  structures.  Because  the  proposed  project  will  involve  the  participation  of 
Professor  Kin  and  will  result  in  specific  and  commercially  available  instrumentation,  I 
am  sure  that  both  the  instrumentation  and  the  theory  and  practice  behind  it  will  be 
incorporated  into  our  Department  of  Engineering's  curriculum.  We  look  forward  to 
Professor  Kin's  active  role  in  the  design,  manufacture,  and  marketing  of  the  proposed 
non-linear  acoustic  instrumentation  for  reliable  damage  assessment. 

Sincerely, 


James 

Chant 


Yackel 

aior 


2200  169th'St  ♦  Hammond  IN  46323-2094  ♦  (219)989-2203  ♦  Fax:(219)989-2581 


3WA  State  University 

r  SCIENCE  AND  TECHNOLOGY 


Geological  Sciences 
253  Science  Hall  I 
Ames,  Iowa  500 11 -32 12 


515294-4477 


terSonix Technologies,  Incorporated 


4615  Post  Oak  Place.  Suite  285 
Houston,  Texas  77027 
phone:  7 1 3-j5£  5-  SBSS,  fax:  7i3-§S3-S£88 
http://ww-v.lasersonix.com  . 

Doctor  Yulian  Kin 
Purdue  University  Calumet 
Engineering  Department 
Hammond,  Indiana  46323 


December  23, 1999 


Dear  Dr.  Kin, 

I  am  writing  this  letter  to  express  our  company’  intention  to  support  development 
of  the  Instrumentation  for  Nonlinear  Acoustic  Nondestructive  Testing.  It  is  my 
understanding  that  Nonlinear  Acoustic  NDE  Technique  opens  a  new  opportunity  in  crack 
and  fatigue  detection. 

LaserSonix  Technologies  is  a  start-up  company,  recently  founded  with  the  aim  to 
develop  a  new  generation  of  imaging  systems  with  sensitivity  and  resolution  that 
supersedes  specifications  of  existing  imaging  modalities.  The  laser  optoacoustic  imaging 
system  (LOIS)  has  two  principal  types  of  applications:  (1)  medical  diagnostic  imaging 
and  (2)  non-destructive  testing  and  evaluation  of  materials  properties:  The  current  goal  of 
LaserSonix  is  to  fabricate  and  market  four  systems  for  the  Opto-Acoustic  Tomography 
(OAT):  (1)  3D  opto-acoustic  tomography  for  breast  cancer  detection,  (2)  2D  endoscopy 
system  for  characterization  of  lesions  in  the  human  GI  tract  from  oral  cavity  to  colon  and 
rectum,  (3)  2D  imaging  system  for  characterization  of  dan  lesions,  (4)  System  for 
characterization  of  properties  of  composite  materials,  ceramics  and  metals  employed  in 
the  aircraft  industry. 

The  sensitivity  and  resolution  of  Nonlinear  Acoustic  NDE  systems  can  be 
dramatically  improved  with  utilization  of  laser  ultrasound  generator  and  wide-band 
acoustic  detectors.  LaserSonix  is  a  world  leader  in  manufacturing  wide-band  acoustic 
transducers  and  also  has  significant  experience  in  laser  generation  of  ultrasound  pulses. 
These  two  system  components  would  be  our  contribution  in  a  collaborative  project  with 
Smart  Material  Design,  Inc. 

We  enthusiastically  support  the  development  of  Nonlinear  Acoustic 
Instrumentation  suggested  by  Smart  Material  Design  and  intend  to  work  together  on  NDE 


Let  it  be  light  that  sounds! 


technique  based  on  laser  excitation  of  probing  ultrasonic  signal  and  applications  of  laser 
shock  as  a  source  of  low  frequency  vibration  signal.  As  a  result  of  our  joint  efforts,  two 
NDE  systems  could  be  developed:  (1)  Nonlinear  Acoustic  Imaging  System  with  laser 
excitation  and,  possibly,  laser  detection,  and  (2)  Laser  Ultrasound  NDE  System  for 
quantitative  material  characterization  at  a  given  location.  LaserSonix  Technologies  will 
assume  participation  in  production  and  marketing  of  the  developed  devices. 

The  following  infrastructure  will  be  employed  by  LaserSonix  Technologies  for  the 
joint  project  with  Smart  Material  Design: 

Compact  nanosecond  YAG-Nd  laser 
Laser  heterodyne  probe.  Model  SH-140 
Front-surface  opto-acoustic  transducer,  model  FST-10 
Tektronix  oscilloscope.  Model  TDS-220  (5  GigaSamples/sec) 

The  cost  of  this  equipment  exceeds  $50k. 

We  estimate  an  annual  labor  cost  for  this  research  and  development  as  $150k  per  year. 


Look  forward  for  our  cooperation. 


Alexander  Oraevsky,  Ph-D. 
CTO 


December  1, 1999 


Dr.  Yulian  B.  Kin 
Smart  Material  Design,  Inc. 

3114  Belwood  Lane 
Glenview,  Illinois  60025 
fex:  847  729  1813 

e-mail-  iMn@calumet.purdiie.edu 

RE*  QUANTITATIVE  NONDESTRUCTIVE  EVALUATION  AND  RELIABILITY  ASSESSMENT 
OF  AGING  AIRCRAFT  COMPONENTS 
Principal  Investigator  Dr.  Alexander  Sutin 

Dear  Dr.  Kin, 

I  am  writing  in  support  of  the  further  development  of  the  above-cited  project,  and  commit  DR S  to  fully 
support  you  in  this  project  I  understand  that  Smart  Material  Design,  Inc.,  Purdue  University,  and 
University  of  Illinois  are  intending  to  extend  their  research  in  nondestructive  characterization  of  materials 
and  aging  structures.  Dynamic  Resonance  Systems,  Inc  is  committed  to  providing  industry  with  the  very 
best  non-destructive  testing  equipment  possible  and  sees  collaboration  as  a  benefit  to  all. 

The  personnel  at  DRS  have  considerable  experience  bringing  new  technologies  to  market-  Bob  Gose,  the 
vice  chairman  of  our  board  of  directors,  is  very  active  in  the  day  to  day  and  long  term  operation  of  DRS 
and  has  40  years  experience  in  the  development  and  utilization  of  advanced  technology.  I  have  been  with 
DRS  since  the  company  was  founded  and  have  been  using  Resonant  Ultrasound  Spectroscopy  for  8  years. 
Data  for  my  PhD  research  at  Colorado  State  University  was  taken  exclusively  with  RUS.  Since  completing 
my  PhD,  I  have  played  an  active  role  in  developing  RUS  as  a  non-destructive  testing  technology.  DRS  has 
been  working  with  Paul  Johnson  at  Los  Alamos  National  Lab  to  develop  non-linear  RUS  as  a  viable 
industrial  tool.  We  would  very  much  like  to  have  DRS  included  in  the  development  of  your  product 

I  would  be  happy  to  discuss  details  with  you  in  person  or  by  telephone. 


Sincerely, 

Frank  Willis,  PhD 
Director  of  Technologies 
DRS 

Fwillis@ndtest.com 


Equipment  /  Personnel  available  at  DRS 

Large  base  of  installed  DRS  Resonance  Ultrasound  Spectroscopy  units  for  both  testing  and  research 
applications 

Microprocessor  development  kits 
Professional  PCB  layout  software 

Visual  C++ programming 

Microprocessor  programming  and  debugging  hardware  /  software 

Established  marketing  organization 

Established  relationship  with  vendors  /  machinists 

Extensive  inventory  of  good  /  defective  manufactured  components 

15  Pentium  /  Pentium  II  computers 

Dedicated  internet  access  /  Web  /  Email  servers 

Oscilloscope 

Function  generator 

Multimeters 

Digital  oscilloscope 

Digital  scale 

Microscope 

Surface  mount  soldering  /  desoldermg  equipment 

Complete  prototyping  capabilities 
Some  machining  capabilities 

MIG  welder  including  aluminum  and  stainless  steel  capabilities 

Large  color  printer  /  plotter 

Color  scanner 
Digital  camera 

StafS^^^Bf1OTgineers,  electronics  technicians  and  marketing  /  adminstrative. 
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